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ABSTRACT 


Middle ultraviolet spectra from 1800 to 3400A are 
analyzed. These spectra were obtained from the March 1990 
rocket flight of the NPS MUSTANG instrument over the altitudes 
meSkm to 315km. The data were compared with computer 
generated synthetic spectra. A least squares fitting 
procedure was developed for this purpose. Each data point was 
weighted using the standard deviation of the means. Synthetic 
Spectra were generated for the following emissions: N, 
wecgard-Kaplan; 9 N, Lyman-Birge-Hopfield; NO vy, 6, and e€; OI 
Pe72A, O11 2470A; and NIT 2143A. Altitude profiles for the 
emissions were obtained. Tentative identification was made of 
the OIII 2853A emission. A comparison of VK and LBH profiles 
demonstrates the process of N, A-state quenching by atomic 


oxygen. 


lil 


TABLE 


I. INTRODUCTION @ a 


A. 


Be 


THESIS OBJECTIVES 


THESIS OUTLINE 


II. BACKGROUND .... . 


A. 


PEL. 


A. 


THE ATMOSPHERE 


OF CONTENTS 


1. General Description .... 


2 SAEMOSDRECL LC AiiaG LOW oe. eee 


ATOMIC AND DIATOMIC SPECTRA 


1. Atomic Bmassions 


2. Molecular Emissions .. . 


DAYGLOW EMISSION FEATURES 2 eeenre 


MOLECULAR NITROGEN 


1. Vegard-Kaplan 


® a ® ® ® 


2. Lyman-Birge-Hopfield .. . 


3. Second Positive 


ATOMIC NITROGEN . 
ATOMIC OXYGEN . . 
NITRIC OXIDE 

i. “Del Ca wee 2. wee 


2. Gamma » & oe 


ney. 


20 


21 


EY 


oir 


“OST kt). Ss ae ae 


Pte EME 3. 6 Sl 


A. 


Je) - 


Cy 


PmiaGauctiOn. “1. . .« . 
Pastreument Deseription .& 
DarcarCollectren J. +. Sa-ae. 


DATA ANALYSIS & RESULTS .. 


A. 


1B) 


ER CHRMMOGUimes . QMS 6 klk lt el 


is 


imercOcuctlOMm +. . >. «oe =. . 
atrttinag Method <9. 5. . . . 


bee oynenetre Spectra... . . 


2. General Programming Technique 
a. Fitting Routine ... . 
b. Weighting Program .... 
PeeWee sos 6 6) #6 @ 4 « «© « 
1. Comparison of Data and Fit . 
a. Nitric Oxide Bands. . 
be Nieregen Bands . ... .. . 
Cc. Oxygen Atomic Emissions 
PIO MOleS as 5 « 5 6 = & « % 
ao erervesOxide... . <a. 2am 
b. Nitrogen Bands... 
Geeetonre EMISSIONS ..<. . . « 
Bee USGUSS POM «  « 2 oe. 6 + © « 


eve OUeENCh@ngG |. 5. . 2 = 


aie 


23 


23 


23 


24 


26 


26 


26 


26 


ai 


30 


30 


33 


34 


34 


34 


34 


38 


38 


38 


43 


47 


48 


48 


be 


Cc. 


Adjustments to Franck-Condon Factors 


Temperature Profile 


VI. CONCLUSIONS 


A. 


B. 


APPENDIX 


APPENDIX 


APPENDIX 


APPENDIX 


APPENDIX 


SUMMARY OF FINDINGS ... . 


TOPICS FOR FURTHER INVESTIGATION. . 


LIST OF REFERENCES 


BIBLIOGRAPHY 


INITIAL DISTRIBUTION LIST. 


Val 


48 


54 


oS) 


oye, 


56 


58 


vs 


96 


105 


110 


216 


218 


220 


LIST OF TABLES 


TABLE PAGE 


2-1 SUMMARY Or SELECTION RULES FOR MOLECULAR 


TRANSITIONS. 5 co, (GSR SIS 6: ee Oe 0s ee 14 
2-2 MOLECULAR ELECTRONIC STATES. ae. S.C. 15 
PemeOmeOmUMnN DENSTIES ~~ 6. « « « « + © «© s « «6 © 8 © « 43 
eee OMG SENT ENSIT TY  VABUES ys SoS eee 6 ee lw we fle 47 
ae ADJUSTMENTS TO FRANCK=CONDON FACTORS .... 4... 51 


18 


LIST OF FIGURES 


FIGURE 

2-1 Temperature profile and corresponding 
atmospheric layers. ...... 

2-2 Atomic processes induced by <* incident 
radiation. o «© « MOC a: << 

2-3 Molecular absorption transitions. 

ZA Graphical representation of molecular 
spectrum. 

FINE Instrument response to molecular transition. 

3-1 N, energy level diagram. 

3-2 Energy level diagram for N’. 

3-3 Energy level diagram for neutral atomic 
oxygen. 

3-4 Energy level diagram for 0”. ........ 

35 Transitwons for nitrile oxide. 

4-1 MUSTANG "hot pixels". 

5-1 "Other Atomics" synthetic spectrum ..... 

5-2 Functions 16) simulate instrument slit 
response.. 

5-3 Nitric oxide synthetic bands at 115km. 

5-4 NO data vs. synthetic spectra at 305, 165, 
105km.. 


Vaan 1 


PAGE 


10 


eZ 


13 


13 


16 


19 


20 


2d 


22 


25 


29 


32 


35 


36 


Comparison of data & LBH synthetic spectrum 
memo L5>5,  lOSKm. 

Viewdatamaxsittt Comparison for 105, 135, 165, 
Lo Skim. 

Data & fit of OI 2972A emission at 115, 155, 
ee okl - 

Data & fit of OII 2470A emissions at 135, 205, 
305km. 

Nitric oxide column density profile. 

Nitrogen LBH intensity profile. 

PMeehaiiiy sprorite For No VK(0,5) Cransition. 
Intensity profiles for OI, OII, OIII, NII. 
Intensity plots for LBH & VK showing quenching 
enke MAK 

Data and best-fitting synthetic spectrum at 
es Kom. 

Data and best-fitting synthetic spectrum at 
305km. 

Unidentified feature at 2955A. 


Temperature profile of analysis. 


1x 


a7 


39 


40 


41 


42 


44 


45 


46 


49 


52 


a 


53 


54 


ACKNOWLEDGEMENT 


This work is incomplete without rendering proper credit to 
my husband, Cleon, whose computer knowledge was constantly 
challenged by my total ineptitude at word processing. His 
refusal to give up on me resulted ina far better composition 


than I could ever have produced alone. 


I. INTRODUCTION 


The composition of the earth’s atmosphere is important 
to many areas of science, such as weather prediction, 
satellite operations, and communications. Properties of the 
upper atmospheric region known as the ionosphere are of 
specific interest to the military, as they strongly influence 
such systems as high frequency (HF) radio communications, 
over-the-horizon radar (OTH), ballistic missile early warning 
(SDI), and the Ground Wave Emergency Network (GWEN). The 
electron density profile, a plot of electron density versus 
altitude for a given time and location, is needed for 
effective use of the ionosphere. Electron density is a factor 
that varies not only with altitude, but with time and 
geographical location as well. In 1986, the Joint Chiefs of 
Staff prioritized the determination of ionospheric electron 
density as fifth on a list of 43 critical geophysical 
parameters for inclusion in a defense satellite (The Joint 
@imets of Staff, 1986). 

Currently, electron density profiles are determined using 
ground-based ionosonde stations. This is very effective 
locally, but is inherently inaccurate for global forecasting 
Since there are only about 20 stations worldwide. A 


satellite-based method of forecasting the ionosphere would be 


globally effective, but to place an ionosonde in orbit is 
unfeasible due to size and power limitations. One alternative 
to ionosondes is to infer electron density profiles from 
spectrographic measurements of the atmospheric airglow. An 
operational spectrograph is small and relatively low-powered; 
it can easily be included in a defense satellite package. The 
Naval Research Lab and the Naval Postgraduate School have a 
joint spectrographic project scheduled to be included in a 
satellite to be launched in 1995. The NRL project is a high 
resolution airglow and aurora spectrograph (HIRAAS) and the 
accompanying NPS effort is the middle ultraviolet spectrograph 
(MUSTANG). The inclusion of the duo represents a significant 
improvement in the global determination of electron density 
profiles over ground-based ionosondes. 

The inference of electron densities from emission 
observations by the HIRAAS/MUSTANG is performed by relating 
measurements of natural ultraviolet emissions to the neutral 
species’ densities. A rocket-launched HIRAAS/MUSTANG 
experiment was first tested in March 1990. The data from the 


launch is the topic of this thesis. 


A. THESIS OBJECTIVES 

This thesis focuses on the analysis of data obtained from 
the March 1990 flight. The analysis extends methods used 
previously on subsets of the data to the entire observed 


range, 1800A to 3400A. All known synthetic models with 


emissions in this range are used in the fitting. The result 
of the fitting algorithm is a set of relative scale factors 
for the various spectra by altitude. The relative scale 
factors are converted into absolute terms and plotted against 
altitude. These absolute plots of intensity and column 
density can be used to infer an electron density profile (see 
eee, Meier ,1991; McCoy,1985). A secondary goal is to 
fetermine the feasibility of an approximation for the N, 
Second Positive spectrum given the unavailability of a 
synthetic spectrum. The determination of the Second Positive 
spectrum would allow an estimation of photoelectron flux, 
Since it is produced by transitions from a state which can 


only be populated by collisions. 


B. THESIS OUTLINE 

The thesis is divided into six chapters. Chapter II gives 
general background information on the descriptions and 
designations of the atmosphere and existing emissions. A 
basic explanation of atomic and molecular spectra is also 
summarized in this’ chapter. Chapter III addresses the 
specific emission features in the wavelength range of 
interest. Chapter IV presents a summary of the experiment and 
the subsequent data reorganization. Chapter V gives a 
complete treatment of analysis of the data and results of the 
study. Chapter VI concludes the paper with a summary of 


findings and suggestions for future study. The Appendices 


contain the computer programs which generated the synthetic 
spectra and the fitting programs, as well as plots of the 


final fits “Geo thnerdata- 


II. BACKGROUND 


A. THE ATMOSPHERE 

1. General Description 

The Earth’s atmosphere is designated by horizontal layers, 
and these divisions are determined by one of two primary 
methods. One method designates the layers based upon 
temperature characteristics; the other uses density to 
categorize the strata. The layers and a typical temperature 


profile are shown in Figure 2-1. 
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Figure 2-1 Temperature profile and corresponding atmospheric 
layers. 


With respect to the temperature designations, the 
troposphere is the region closest to the Earth’s surface where 
temperatures decrease with altitude, typically from about 300K 
at the surface to 220K at 10 kn. The point at which 
temperature begins to rise again is called the tropopause and 
it indicates the start of the next region, the stratosphere. 
The stratosphere is a relatively warm layer due to absorption 
of solar radiation(primarily 2000-3000A) by ozone. The top of 
the stratosphere 1s marked by the stratopause. At this point 
the mesosphere begins. The temperature again begins to 
decrease with altitude until a low of about 100K occurs at 
around 85 kilometers. This is the mesopause, which marks the 
beginning of the next layer, the thermosphere. In the 
thermosphere, temperature rises steadily with increasing 
altitude until reaching an equilibrium temperature between 
1000K and 2000K, depending on solar activity. 

Studying the atmosphere in terms of density gives rise to 
a different set of atmospheric layers. Under this convention, 
the lowest 80 kilometer layer of the atmosphere is called the 
homosphere. In this region, vertical mixing of the air 
currents provides a uniform relative density ratio between 
various species. The primary components are nitrogen, oxygen, 
and argon in a ratio of 78:21:0.9. The remaining 0.1 percent 
is primarily helium, hydrogen, and carbon dioxide. The 20 km 
from 80 to 100 is designated the barosphere, or tropopause, 


and it is a region of transition above which very little 


vertical mixing occurs. Above the barosphere, from 100 km to 
around 1000 km, lies the heterosphere. Due to the absence of 
mixing in this region, the ratio of constituents varies with 
altitude. At extremely high altitudes of 1000 kilometers or 
so, the atmosphere becomes so tenuous that particle 
trajectories are essentially parabolic orbits. This "outer 
atmosphere" is called the exosphere. 

The data analysis in thiS paper is concentrated in the 
lower heterosphere, the first 200 km of the thermosphere, 
where there 1S an obvious positive temperature gradient and 
the convective mixing iS minimal. This region is also termed 
the ionosphere. 

The ionosphere is so termed due to the ionization of the 
atmosphere by absorbed solar radiation. The actual height and 
thickness of the four layers(D,E,F1,F2) vary with solar 
activity and atmospheric composition. The D region is 
typically from 70 km to 90 km during the day; the E region is 
typically from 90 to 130 km; Fl generally exists from 130 to 
160 km during the day; and F2 is above 160 km. Each of these 
regions is formed by ionization due to solar radiation at 
different wavelengths. For example, longer wavelength 
radiation can penetrate deeper into the atmosphere. Since 
solar activity changes daily, seasonally, and yearly, the 
intensity of sunlight at a particular wavelength and hence the 
depth of penetration is variable. Ionization depends on the 


composition of the neutral atmosphere as well. Combining 


these two factors of composition and solar activity, one can 
see that the flux of radiation is greatest at high altitude 
but density is lowest and vice versa. The maximum ionization 
occurs somewhere between these two extremes, and the result is 
an overall layer known as a Chapman layer. The composition 
and formation of the ionosphere is quite complex and is not 
unique to Earth. An excellent discussion of planetary 
1onospheres is written by Chamberlain (1978). 

2. Atmospheric Airglow 

Airglow is the naturally occurring radiation in the 
Earth’s atmosphere. The daytime emissions, called dayglow, 
are primarily driven by photoexcitation. Also contributing to 
dayglow are emissions which result from photoelectron impacts 
and, to a lesser degree, photochemical reactions. During 
daylight hours photochemical reactions are negligible compared 
to the abundance of emissions due to sunlight, but the same 
reactions are significant at night since there are obviously 
few direct solar photons to stimulate transitions. Because the 
data analyzed in this work were obtained during the day, 
photochemical reactions are not specifically addressed in this 
work. 

The primary source of dayglow emissions is 
photoexcitation. Photoexcitation is the absorption of a 
photon as a result of inelastic collision between a photon and 


an atom or molecule to raise the particle to an excited state. 


When the photon collides with the particle, the photon is 
absorbed and the particle jumps to a higher energy level. 
When the particle drops back to a lower energy state, it emits 
a photon with energy equal to the transition. This process 
will be further described in the next section on atomic and 
molecular emissions. 

Although not as important to dayglow as photoexcitation, 
photoelectron collisions are significant because they can 
produce transitions which are otherwise forbidden. 
Observations of forbidden transitions can be used to find the 
photoelectron flux, since only photoelectron collisions could 
excite the particle to the upper state. Given the importance 
of whether a transition is allowed or forbidden to 
calculations of photoelectron flux, a knowledge of basic 
selection rules is desirable for the analysis of dayglow 
emissions. 

It is also worthy to note that some of the forbidden 
transitions which are prominent in the low density unconfined 
gas of the thermosphere cannot be produced in earthbound lab 
experiments (Rees,1989). Such emissions, seen in dayglow data 
collected by space-borne equipment, do not match any of the 
spectra which have been observed using lab or ground-collected 
data. These emissions are identified with transitions of 
Known elements, by wavelength matching. One needs to be 
familiar with spectral nomenclature and selection rules in 


order to understand the explanations of these identifications. 


B. ATOMIC AND DIATOMIC SPECTRA 

1. Atomic Emissions 

The primary atom-photon interactions are shown in Figure 
2-2. In elastic scattering, no energy transfer takes place 
between the photon and the atom during collision, although the 
photon may change directions. In inelastic scattering, the 
incident photon may not only change direction, but may impart 
energy to the atom. In this case, the atom is excited to a 
level AE above its original state and the deflected photon 
will have an energy AE less than its incident energy. A 
photon with energy AE will be emitted when the atom relaxes to 
the initial state. Resonance scattering occurs when the 


incident photon is absorbed, raising the atom to an excited 


OG 


a) Elastic scattering b) Inelastic scattering c) Resonance radiation 


=F . = 
fy “a, 
Se 


d) Fluorescence e) Collisional De-excitation 





Figure 2-2 Atomic processes induced by incident radiation 
(Brehm & Mullin, 1989). 
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state. A photon with wavelength equal to the incident photon 
is emitted when the atom falls back to its initial state. 
Fluorescence is the phenomenon which occurs when the excited 
atom undergoes several downward transitions in returning to 
the initial state, giving a cascade of photons. An atom can 
also return to the ground state by colliding with another 


particle, thus returning to the ground state without emitting 


a photon. This 1s called collisional de-excitation or 
quenching. 

Similar ee) photon-atom collisions, electron-atom 
collisions may be elastic or inelastic. In an electron-atom 


collision, the scattered electron loses kinetic energy, rather 
than actually being absorbed, when the atom is excited to a 
higher energy level. Photon emissions occur when the atom 
relaxes to the ground state. 

2. Molecular Emissions 

Molecular spectra may occur as the result of electronic 
state changes, but they may also result from purely 
vibrational and rotational state changes. Each of these types 
of transitions has typical energy gaps, with electronic being 
on the order of 10 eV, vibrational of .1 eV, and rotational of 
-001 eV. These three types of transitions are illustrated in 
Figure 2-3. As seen in the figure, for a given electronic 
energy level there are many vibrational energy levels. 


Similarly, for a given change in vibrational energy there are 


ude 


Energy 





Figure 2-3 Molecular absorption transitions: 
Transition I is pure rotational 
Transition II, vibrational change Ap=1 
Transition III, electronic change of state 
(Brehm & Mullin, 1989). 


many possible initial and final rotational levels. AS a 
result of the many vibrational and rotational levels which can 
be populated for any given electronic transition, the emission 
feature is composed of many closely spaced lines. Figure 2-4 
shows a molecular spectrum, exhibiting a distribution of 
emissions. When the separation of component lines is smaller 
than instrument resolution, the transition may appear as a 
hump, as shown in Figure 2.5. This is the case for the data 


studied in this paper. 


2 


The study of quantum 


mechanics explores 
transitional possibilities, 
and develops a set of 


guidelines for those emissions 


one expects to observe These 





so-called selection rules 

require that rotational 
Figure 2G Graphical 

transitions only occur for representation of molecular 
spectrum. 


4j=t1, where j is the total 
feud ah t um mn wm Se -r 
m@—0,1,2,3;,-+-)- Vibrational 
transitions occur within the 


Gonrines of “selection “rules 


aj=t1 and av=tl1, where v is Figure 2-5 Instrument response 
to molecular transition. 


the vibrational quantum 
number (v=0,1,...j}). These rules are not strictly adhered to, 


as they are based on electronic dipole behavior, and some 


molecules exhibit quadrupole and magnetic dipole 
characteristics. There are also cases where, though only 
slightly probable, emissions occur when aAv>ti. Electronic 


transition selection rules are complex, and a full treatment 
1s given to the development and explanation of them by Green 


and Wyatt (1965). A summary of these is shown in Table 2-1. 
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TABLE 2-1 SUMMARY OF SELECTION RULES FOR MOLECULAR 
TRANSITIONS (Green & Wyatt, 1965). 


Electric Quadrupole (E2) 


Total inversion + <3) 4 
Electronic 
inversion g°9 ON i ed 
0*,+1 (0 + 0) QO, +1, +2} 
0 
0 


Geri, £2 


El ee 
eel +9 
0, £1§ +2 





*NotO—O fdt0-0,0 -1,4-} 
TIf0— OthenAy=0. §Not for 5 — 3 transitions. 


i 
Derivation of these rules is not essential to this data 
analysis, but an understanding of the corresponding molecular 
nomenclature is necessary in order to understand the 
discussion. 

The component of total angular momentum which lies along 
the internuclear axis of a molecule is designated A, and is 
used to label molecular electronic state. Each state is given 
a capital Greek letter depending on its value of A. These 
designations are shown in Table 2-2. The spin of the molecule 
is designated = (which should not be confused with the 2& 
symbol for the A=0 state). It can take on values from -S to S, 


where S is the resultant of individual electron spins in the 
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TABLE ie MOLECULAR 
ELECTRONIC STATES. 
molecule. Mae 1 Cr ry , WCB FUE 


ear ils is written as a left DESIGNATIONS OF 
MOLECULAR STATES 
superscript on the state letter. 
Greek Ok 
The ground state of a diatomic am a 
molecule is labelled with a 


leading xX. Successive excited 


states of the same multiplicity 





are labelled A,B,C, and so forth. 
Excited states of different 
multiplicity are alphabetically 
designated by small letters. If these labels are strictly 
followed, then in general the selection rules’ allow 
transitions between like-labelled states (eG. 
capitalScapital). As with any rule, there are exceptions. 
For instance, N, diagrams have capital and small letters 


assigned conversely. 
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III. DAYGLOW EMISSION FEATURES 


A. MOLECULAR NITROGEN 

Figure 3-1 is an N, energy level diagram including the 
three transition bands of interest in the 1800-3400A range. 
The states are indicated by 
molecular notation, and the 
transitions are identified by 
name. 

1. Vegard-Kaplan 


The Vegard-Kaplan system 


Second 
positive 


is so designated because 


~ 
2 
o 
Ww 
> 
0 
de 
y 
= 
UJ 


Vegard discovered the system 
in 1932, and Kaplan determined 
the specific wavelengths 
involved a few years later. 
The emissions range from 1250A 


to 5325A and represent the 





3 yer 1 + : . ; 
APL, 7X 2Z, transition. In this Figure 3-1 N, energy level 


Eransa Elen eae total SH uchepecc un (iSdIl a ENS) « 


transformation is ++, the electronic inversion is u-g, Ar=-1, 
and AA=0. From the molecular selection rules summarized in 


Table 2-1, this is a forbidden transition, yet it has been 
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observed as an electric dipole transition. het horough 
discussion of the band system is given in Meier (1990). 

A-state nitrogen molecules are produced by cascading 
downward from the B and C stateS or as a result of 
photoelectron collisions. Because the relationship of 
emission bands to photoelectron flux is complicated by 
cascading N, molecules and quenching effects, the V-K band is 
not generally used to infer photoelectron flux. 

2. Lyman-Birge-Hopfield 

The Lyman-Birge~-Hopfield system is the name given to the 
a'TI,>x's,* transition of N,. The emission bands arising from 
this transition extend from 1160A to 3020A and are due to both 
electric quadrupole and magnetic dipole transitions. Unlike 
the V-K system, the upper state of the LBH system can only be 
populated by photoelectron collision excitation. This system 
can therefore be used to directly infer the photoelectron 
fe eux « The LBH bands found in the data were particularly 
significant in the wavelength region from 1800A to 2400A. An 
intensive investigation of the LBH system as it pertains to 
the MUSTANG data was performed by Mack (1991). 

3. Second Positive 

The Nitrogen Second Positive system (C°II>B‘II,) extends 
from 2680A to 5460A. It is a particularly important 
transition because the upper level is not directly populated 


PiLoughepmotocxeitation, but only by collision. Thus a 
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measurement of the Second Positive intensity profile can be 
used to infer the photoelectron flux. This result can be 
compared with those obtained from the  LBH bands. 
Unfortunately, there is no synthetic spectrum model developed 
for the Second Positive. It was originally intended to 
determine the N, Second Positive spectrum by subtracting known 
contributions from the data and assuming the remaining data 
was Second Positive. This means of estimating a spectrum was 
found to be impractical due to inaccuracies in the data caused 


by data drops. 


B. ATOMIC NITROGEN 

Figure 3-2 shows the transitions which have been observed 
in the airglow due to N’. Two transitions are in the 
wavelength range of this study. The 3070A emission falls in 
the same range of wavelengths as the N, Second Positive 
spectrum, so it could not be definitively seen in the data. 
The emission shown at 2143A was considered in the synthetic 
fit to the MUSTANG data. It is actually a doublet at 2143.55A 
and 2139.68A with a ratio of 100:58, respectively. This 
intensity ratio was determined by Bucsela and Sharp (1989) 
from dayglow observation. Theoretical and laboratory 
estimates of relative line strength range from 100:60 to 
100:23. Since the data source for the Bucsela and Sharp ratio 
waS Similar to the experimental data being analyzed, the 


100:58 relation was used in the synthetic model. 
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Figure 3-2 Energy level diagram for N° (Rees, 1989). 


C. ATOMIC OXYGEN 

Figures 3-3 and 3-4 diagram the atomic oxygen emissions 
which have been viewed in the airglow. Two of these emissions 
fall in the wavelength range explored by this analysis. 

The neutral oxygen atom emits a photon with wavelength 
2972A when it falls from the 'S state to the °P ground state. 
This transition is depicted in Figure 3-3. Singly ionized 
oxygen, shown in Figure 3-4, produces a line at 2470A when it 


eels ion the -P state fo the “Ss ground state. 
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Figure 3-3 Energy level diagram for neutral atomic oxygen 
(Rees, 1989). 


D. NITRIC OXIDE 

Nitric oxide emissions, shown in Figure 3-5, dominate much 
of the dayglow data. There are four primary bands in the 
range 1800A - 3400A, vy, 6, B, and €. The 6 bands are several 
orders of magnitude less intense than the other emissions, so 
they have been neglected. 

1. Delta 

The C*II-x‘Il transition band system, designated 6, ranges 
from 1621A to 4288A (Barth,1965). Using standard notation, 
the upper electronic state’s vibrational level is designated 
v’, and the lower state vibrational level, v". The 6 


transition band includes vibrational levels corresponding to 
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Oe re eemande 0, 1,...23. Nitric oxidewwill dissociate 
Be v'>O for Cll, so v’ is limited to v’=0 and v"<14 in 
creating the synthetic 6 spectrum. 

2. Gamma 

The y bands of nitric oxide extend from 1873A to 6126A. 
The emissions arise from the A*=*-x*°II transition, and the 
ground state has v"=0,...23 as previously stated. The 
dissociation energy for nitric oxide is lower than the energy 
of the v’=4 level in the A state, so only transitions arising 


from v’=0,1,2,3 contribute to the synthetic y spectrum. 
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Figure 3-4 Energy level diagram for O° (Rees, 1989). 
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3. Epsilon 

The nitric oxide € system 
is the result of molecules 
transitioning from D*5*=x°II. 
The wavelengths of emissions 
range from 1800A to 2385A. D- | Ogawa 2 


state molecules populate 


Vibrational levels v’=0,1,2. 





Figure 3-5 Transitions for 
nitric oxide (McEwan, 1975). 
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IV. THE EXPERIMENT 


A. Introduction 

The middle ultraviolet spectra analyzed in this thesis 
were obtained from the NPS MUSTANG spectrograph flown on a 
NASA sounding rocket launched March 30, 1990, from White Sands 
Missile Range, New Mexico. The spectrograph was launched on 
a Terrier boosted Black Brant, which carried the payload to an 
apogee altitude of 320 km. The spectrograph was looking in 
the anti-solar direction with an observation zenith angle of 
90°. From 320 to 180 km on the downleg no useful data were 


obtained due to the payload attitude. 


B. Instrument Description 

The MUSTANG (Middle Ultraviolet Spectrograph) is a 
modified 1/8 m Ebert-Fastie spectrograph with a 1/8 m 
telescope. The entrance slit is 5mm x 140 pum and a 10A 
wavelength resolution is provided by a 1200 g/mm grating 
ruling. The detector assembly consists of an ultraviolet-to- 
Visible image intensifier using a quartz input window, a CsTe 
photocathode,and a fiber optic output window. The detector is 
an array of 512 photodiodes giving voltage outputs 
proportional to intensity using a 50 msec integration period. 
This voltage output is fed serially to an A/D converter and 


then to rocket telemetry via a FIFO shift register. 
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Calibration/sensitivity curves were obtained using a 
deuterium lamp for the wavelength range 1800-2800 A and a FEL 
tungsten lamp for 2400-3400A. A more precise calibration of 
the 1800-2000A range, which is affected by 0, absorption, was 


performed in a vacuum chamber (Mack,1991). 


C. Data Collection 

Over 8000 spectra were collected by MUSTANG. These 
spectra were then averaged into 25 altitude bins, disregarding 
spectra which were contaminated by the  instrument’s 
electronics. Each of the 25 resultant spectra covers a 10 km 
layer and is identified by its central altitude. 

Because of a timing difference between clocking and 
reading functions, only 480 of the 512 photodiodes’ responses 
were recorded. It was first determined that every seventeenth 
word waS dropped, and these blanks were inserted aS averages 
of the two adjoining data values (Clayton,1990). This still 
did not produce a suitable wavelength match with known 
emissions, so a mathematical spline was performed on the data 
in order to match a uniform wavelength grid (Mack,1991). 
Although this provided an excellent fit to assumed models, it 
had no physical justification. Laboratory testing was 
performed by Quint (1991) to determine which values were lost 
in the timing error. His results accounted for the 32 
expected data dropouts, but two additional dropouts remained 


unexplained. These two shifts were inserted at positions 
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which were consistent with the findings of Quint and also 
provided the best agreement with three previous studies. It 
should be pointed out that this was not confirmed with the 
electronics. The corrected spectra are delineated in Appendix 
A, noting each dropped pixel. At these spots an average of 
the adjacent values was inserted into the original 480 data 
points, to obtain an array of 512 values. 

The adjusted spectra were then divided by instrument 
Sensitivity. The first six 
values were set to zero, since 
the uncertainty in sensitivity 
of the instrument at the 
endpoint is extremely large. 
It was found that one further 
correction had to be made to 
the data. In both pre-launch 
and post-launch high voltage 
dark count spectra, two 
voltage spikes existed. These 


hot spots occurred at points 





corresponding to 2850A anda 


Figure 4-1 MUSTANG “hot 


3373A, as shown in Figure 4-1 Pixels" plotted as intensity 
(in R/A) versus wavelength. 


These dark count hot spots 
were converted to R/A and then subtracted from the averaged 


spectra before analysis. 
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V. DATA ANALYSIS & RESULTS 


A. TECHNIQUE 

1. Introduction 

a. Fitting Method 

Synthetic spectra were fit to the data using a least 
squares fit process (hereafter referred to as GRIDLS). It was 
developed by Bevington (1969) and applied in limited portions 
of the wavelength range by Clayton (1990), Andersen (1990), 
and Mack (1991). This method is an extension of a basic 
least-squares regression, minimizing error between the fit and 
the data as measured by the variable y* with weighting for 
confidence in the data. The program fits a linear combination 
of synthetic spectra to the data, with each synthetic spectrum 
having an independent coefficient, or scale factor. The 
program is adaptable for any number of desired parameters. 
The final analysis used nine parameters, which will be 
discussed in the Results’ section. The fitting progsam 
(GRIDLS) and its accompanying subprograms are in Appendix B. 

The output of the program is the scale factors for the 
linear combination of models. What each models’ scale factors 
represent depends on the normalization dimensions used in 
creating the corresponding synthesized spectrum. The model 


spectra do not all have the same dimensional units. In the 
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case of atomic spectra, the scale factors are absolute 
intensity measured in Rayleighs. For the VK and LBH bands the 
scale factors are a measure of relative intensity and must be 
converted to absolute values to obtain an intensity profile. 
The Nitric Oxide scale factors are column density, in 
molecules per square centimeter. 

b. Synthetic Spectra 

Previous analysis of the March 1990 MUSTANG data includes: 
an analysis of NO at 2000-2500A (Clayton,1990) ; analysis of OI 
2972A and OII 2470A emissions (Andersen,1990); and a study of 
LBH and NO at 1800-2100A (Mack,1991). An analysis combining 
all the contributing spectra and spanning the full MUSTANG 
wavelength range had not been performed, but previous work 
provided many of the computer programs as well as "initial 
guesses" for the fitting routine. Other models and programs 
arose as a result of the fitting analysis. Computer models 
written by Cleary (1986) were used for the generation of 
Nitrogen Lyman-Birge-Hopfield and Nitric Oxide y, 6, and € 
band emissions. The Nitrogen Vegard-Kaplan model was written 
by Siskind and Barth (1987). 

The synthetic spectra for individual atomic emissions were 
created by convolving appropriate impulse functions with the 
slit function of the instrument. As already discussed, in the 
case of the N* 2143A doublet a pair of impulse functions in 


the ratio 100:58 was used. The previously identified atomic 
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contributions to the MUSTANG spectrum include OI 2972A, OII 
2470A, and NII 2143A. 

One other atomic spectrum, a prominent line emission at 
2853A at higher altitudes, has been included in the fitting 
routine as a result of observation. Elements with emissions 
at that wavelength include nickel, copper, manganese, iron, 
sodium, cobalt, and oxygen. Of these, oxygen is the most 
abundant in the atmosphere where the emission is seen, above 
200km. Doubly ionized oxygen has an emission at 2853A that 
although not previously identified in dayglow data, is 
classified "predicted" by Kelly (1979). 

A collection of atomic emissions observed during the 
fitting process, but not verified by other literature, was 
classified as a separate spectrum. As discussed in Chapter 
II, an atomic emission is a Single line vice the "hump" of a 
molecular transition. Even after convolution with the 
instrument slit function, the form of an atomic emission 
feature is distinctive from a molecular band. There were 
several areas in the data which matched the form at a 
sufficient number of altitudes as to suggest an atomic 
emission. These wavelengths were then checked against tables 
for possible origins and included in the synthetic spectrum if 
a likely emission existed, based on atmospheric composition 
and transition probability. The resultant synthetic spectrum 
(designated "Other Atomics") is shown in Figure 5-1. The 


program used to create this spectrum, along with possible 
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"Other Atomics" synthetic spectrum 


Figure 5-1 
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identifications for each emission, iS in Appendix B. None of 
these emissions occurred with sufficient intensity and 
regularity to merit being considered as a Separate parameter. 
They are included only to facilitate the fitting procedure. 

2. General Programming Technique 

a. Fitting Routine 

The GRIDLS program which succesatullly aieecnedis points of 
data with three parameters (Mack, 1991) was overloaded by the 
inclusion of all 512 points of wavelength and nine parameters 
(i.e., the scale factors for V-K, OI, OII, OIII, LBH, NO, NII, 
Other Atomics, and Background). Several modifications were 
made to improve the data handling capabilities of the progran, 
and other characteristics were accepted as limitations. 

Previously, the program would allow negative scale 
factors, effectively allowing subtraction of a spectrum. It 
also would enter an infinite loop if the scale factor equalled 
zero. These problems were eliminated in the program by 
requiring that zero be the minimum scale factor allowed, and 
by performing a Simple global minimum check to prevent a zero 
scale factor from causing an infinite loop. 

The program had great difficulty converging to a value if 
the "initial guess" for the scale factor was far from the 
minimum deviation value. This problem was not corrected, but 


a test-plotting program was written to allow one to visually 
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match the synthetic curve to the data so as to obtain 
reasonable first guess values. 

The standard function of instrument response was an 11- 
element symmetric slit function (Clayton, 1990; Andersen, 
m90; Mack, 1991). When convolved with a 512-element 
synthetic spectrum, this slit function would sufficiently 
represent the instrument’s response to the modeled molecular 
spectrum. However, in the case of an atomic emission, which 
occurs as a line at a particular wavelength vice a band of 
lines, the slit function could shift the apparent peak as 
much as 1.7A. The symmetric function centered the response on 
the applied pixel, a 3.133A range, instead of at the 
particular desired wavelength. A new set of five slit 
functions, shown in Figure 5-2, was developed to correctly 
model atomic emissions, allowing for emissions which are not 
centered on a pixel. 

A major limitation of the previous version of the GRIDLS 
fitting routine was that there existed no absolute measurement 
of "goodness of fit." Although the program minimized a x° 
value, the values were not weighted and thus could not be 
compared between runs or altitudes to determine goodness of 
fit as parameters were tweaked. Creating a proper weighting 
for the data evaluation so that these comparisons could be 


made waS a major task. 
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Figure 5-2 Functions to simulate instrument slit response. 
Delta function is placed: a)centered b)% pixel right of 


center c)% pixel right d)% pixel left e)% pixel left. 
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b. Weighting Program 

The GRIDLS program allows for statistical weighting of the 
data, instrumental weighting, or no weighting. Instrumental 
weighting was chosen for its advantage of stressing datapoints 
which have small uncertainties. The first step in producing 
instrumental data weights was to determine the standard 
deviations of the mean (sdom) for every data point in each 
altitude bin. A second source of instrumental uncertainty was 
instrument response. It was determined that the detector had 
a linear response for most output values, but reached 
Saturation above 900 decimal units (divided by instrument 
sensitivity to get R/A), so data points in the saturation 
range were given weights based on the sdom being equal to the 
data value itself. The exception to this was when the data 
value was so small as to give a large weight to the point 
(since weight «sdom'), in which case the weight was assigned 
1/10 the largest weight factor for that altitude (Cleary, 
private communication, 1991). Additionally, the data points 
which were merely averages of adjoining values due to lost 
data points were given weights based on the sdom being equal 
to the data value. Since the sdom and saturation points are 
different for each altitude, a set of values for weighting of 
the 512 data points was stored for each data spectrum. 

The data pixels corresponding to nitrogen Second Positive 
emissions had no synthetic input being matched to them. The 


pixel corresponding to each wavelength of these emissions as 
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determined by Barth (1965), plus a buffer of +1 pixel to 
account for instrument response, were given essentially zero 
weight. The final weight program, designated WEIGHTS2, is 


Appendix C. 


B. RESULTS 
1. Comparison of Data and Fit 
a. Nitric Oxide Bands 

In the plots of the fits, nitric oxide is viewed as one 
entity rather than plotting the y, 6, and e€ bands separately. 
This is for ease in reading the plots. Figure 5-3 shows the 
components of the nitric oxide synthetic spectrum for 115km. 
The bands are not equal contributors, as the & band has a 25% 
fluorescent efficiency factor. 

Figure 5-4 1s a plot of the observed data and its 
corresponding fit. The data is indicated by the solid line, 
and the fit obtained from the computer procedure by a dotted 
line. The dashed line is that portion of the synthetic fit 
which is nitric oxide. From the plots in Figure 5-4, one can 
see the importance of nitric oxide with decreasing altitude. 
In the 105km plot, nitric oxide accounts for almost all of the 
dayglow from 1800-2600A, but at 295km the NO is less 
prevalent. 

b. Nitrogen Bands 
Figure 5-5 shows an intensity plot (in R/A) of the LBH 


contribution to the fit, along with the observed spectrum, for 
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Nitric oxide synthetic bands at 115km. 


Figure 5-3 
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three different altitudes. From these plots, one can see that 
the Lyman-Birge-Hopfield band of N, emissions dominates the 
wavelength range 1800-2000A. 

Figure 5-6 shows a Similar altitude progression for the 
Vegard-Kaplan system, emphasizing the (0,5) transition at 
2605A. VK emissions are prominent from 2200-3300A, and in 
general all the emission features show the same characteristic 
growth as exhibited by the (0,5) transition. 

c. Oxygen Atomic Emissions 

The OI 2972A emission as fit to the data at three 
different altitudes is shown in Figure 5-7. The poor fit of 
this feature is addressed in the Discussion section. The OII 
2470A line, shown in Figure 5-8, becomes more important with 
increasing altitude. 

2. Profiles 

a. Nitric Oxide 

The column density profile for Nitric Oxide is shown in 
Figure 5-9, with the corresponding values given in Table 5-1. 
Values obtained from upleg data are indicated with squares and 
Gdownleg values by triangles. There is a systematic 
uncertainty Bef ils The values obtained show general 
agreement in both magnitude and slope with the column 
densities obtained by Mack (1991), Clayton (1990) and Andersen 
(1990), all of whom analyzed subsets of the MUSTANG data. The 


maximum and minimum values best agreed with those obtained by 


38 


Data & Fit at 105km (downleg) Dota & Fit ot 135km (downleo 


Intensity(R/A) 
Intansity(R SA) 





ene 


9700 7800 9900 2600 2700 2600 2900 









Wavelength(A) Wovelength(A) 
Data & Fit ot 155km (uplea) Dato & Fit at 195km (uplec 
200 300 
—— Data 
—— Fit — Data 
400 : - --Vk = aa | 
500 


—— 

mn 

=") 
—— 


Iintensity(R/A) 
Intensity(R/A) 
BS 
=) 
> 
fA ia 








20 an 
eo: 
| pS 
OL: 
2600 2700 2600 2300 2600 2700 2800 2900 
Wavelencth(A) Wavelength(A) 


eritttitriias lisa 
; (a a 


= F 


Pcuee 5-6 “Vk data & fit Comparison for 105, 135, 165, 
195km. 


oo 


Intensity(R/A) 


Intensity(R/A) 


Figure 5-7 
275K 


Intensity(R/A) 


800 


600 


200 


300 


200 


40 


N) 
O 


8) 
2940 


Data. Soffit -of “Ol 2972ANemissilon. ae 


Dato & Fit at 115km (downleg) 


Dato 
“as Fit 
=) Ol 29724 


2960 2980 
Wavelength(A) 


Data & Fit at 155km (upleg) 


ae eee Data 
aes a | 
So ee ee Ol 2972A i 





2960 2980 
Wovelength(A) 


Data & Fit at 295km 


2960 2980 
Wavelength(A) 


40 


5000 


 ERERTPIETPIPEPTIPESLOSINISIDTEI RTT INEY 


300 






155, 


Datauccsen sat 15okm (downlec 


Intensity(R/A) 





2400 2450 2500 Zoo0 
Wavelencth(A) 


Data & Fit at 205km (upleg) 


S00 
Se... Data 
2907 ie . Fit at 
aa O! 2470A 
= 209 
ieee 
oS 
= 150 
w \ 
D a 
= 400 see: | 
aan : 
| 
SS = 
r | 
| L 
9 t= L 
2400 2450 2500 2a5C 
Wavelength(A) 
Data & Fit at 305km (uplea 
130 
Sa Dota 
— ai 
ee = OH 2 470A | - 
oT 
i of 
[8 ‘ 
es | 
"A | 
LS vie 
S : 
= ise aaa. 
| 
! 
ne L 
O ! 
2400 2450 ZIC0 Zo00 


Wavelength(A) 


Figure 5-8 Data & fit of OII 2470A emission at 135, 205, 
305km. 
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Figure 5-9 


42 


Mack (ano Sas), with the 


intermediate values best 


matching the fit obtained by 
Clayton (1990). The 
Significance of this agreement 
is that mL AL 


a spectrum 


feel te eer Ch procedure (li.e., 
including all species and all 
wavelengths) is not required 
to accurately determine the 
Nitric Oxide density. 
b. Nitrogen Bands 

The scale factors for LBH 
bands were converted to 
absolute intensity using the 
algorithm and computer 
programs in Appendix D. 


is shown in Figure 5-10, 


intensity as it includes all a-X transitions. 


TABLE 5-1 NO COLUMN DENSITIES 


| ALTITUDECKm 
) 
*downleg 





COLUMN DENSITY 
(molecules/cm’) 


1.908e+16 
1.681e+16 
1.186e+16 
8.785e+15 
5 .883e+15 
2.427e+15 
3.596e+15 
1.938e+15 
2.299e+15 
1.779e+15 
1.520e+15 
1.332e+15 
9.850e+14 
8.222e+14 
6.321e+14 
5.062e+14 
4.003e+14 
3.248e+14 
2.554e+14 
2.103e+14 
1.665e+14_ 





1.158e+14 
9.824e+13 
6.240e+13 








The resultant LBH intensity profile 


and 1S a measurement of total LBH 


The intensity 


profile for VK can only be calculated for a particular v’-v" 


band, 


and Figure 5-11 shows this calculation for one of the 


more prominent systems, the (0,5). 
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Figure 5-10 Nitrogen LBH intensity profile. 
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Figure 5-12 Intensity profiles for OI, OII, OIII, NII. 
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c. Atomic Emissions 


The intensity profiles for OI 2972A, OII 2470A, OIII 


2853A, and NII 2143A 
are shown in Figure 
5-12, with the values 
maven in Table 5-2’. 
The OI 2972A profile 
shows a distinctive 
peak at about 150km; 
mnemO’ =2470A profile 
shows a broad peak at 
about 250km; the OIII 
profile suggests a 
peak at around 300km; 
and the NII reaches 
its maximum around 
BS Okm. 

The OI 2972A line 


emission profile was 


TABLE 5-2 ATOMIC INTENSITY VALUES 


| 105" | 1878.70 | 7816.57] 0.00] 934.11 
ise eee eso: oniadey> | OsCON me Misr95" 
POSE eae ee 
| 135 | 1187.29 | 1660.38] 0.00 | 3106.81 
| 145* {| 1652.85] 1633.97] 0.00 | 2824.51 | 
| 155 | 1777.33 | 2065.98] 0.00} 1824.37 
| 155" | 1828.75 | 1323.66] 0.00 | 3469.86. 
| 165 | 1623.30] 1671.16} 0.00] 1715.49 | 
_ 165 | 1692.26] 1264.87] 0.00] 2559. 10 | 
| 175 | 1351.22] 1436.91] 0.00] 1431.54 | 
Pp 175*{ 1448.04] 1127.43 | 0.00] 2175.71 | 


| 185 | 1045.54] 1478.36 | 0.00] 1243.45 
| 195 (| 907.09] 1535.49] 0.00] 1238.92. 
| 205 | 677.06] 1517.33 | 0.00] 865.34 
| ats | 461.95] 1492.67] 0.00] 786.54 
| 225 | 404.27] 1523.80] 0.00] 654.12 | 
a5 | ese is 25955i | jonool|| 563259) 


| 245 | 353.83] 1540.26| 0.00] 496.13 
255 | 240.10} 1522.03] 0.00] 436.25 
348.22 
PE 

| 169.77] 1383.26 | 25.48 | 
295 | 156.30| 1350.18] 30.51 | 
| 147.63 | 1342.65 | 32.31 | 





found to have very little agreement with previous analysis. 


The maximum found by this analysis was 50% of that found when 


fitting only the OI 2972A line (Andersen, 1990). 


This is due 


to the uncertainty of the wavelength for the data, since two 


of the data dropouts were in this region. 


The overall 
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uncertainty in the OI profile should reflect the data drops, 
thus caution is urged for any subsequent analysis of this 
emission profile. 
3. Discussion 
a. V-K Quenching 
Figure 5-13 shows an overlay of the nitrogen LBH and VK 
profiles, with the VK scaled up by a factor of 9.5 for 
comparison. The most interesting aspect of the nitrogen 
profiles is that they have roughly the same slope above 200 
kilometers, but the VK intensity drops off rapidly below this 
height. This is due to the quenching effect of N5, NO; Meme, 
and N at the lower altitudes. From auroral studies, atomic 
oxygen is thought to be the primary quenching agent with 
molecular oxygen playing a role below 135km (Sharp, 1971). As 
previously discussed, the upper state for the VK emissions is 
the N, A-state. The A-state to ground-state energy for N, 
corresponds to an excitation energy for atomic oxygen. Below 
200 km, the O density is high enough that during their 
metastable Jufetimemin the excted@state, N molecules comme 
with and excite oxygen, thus losing their energy without 
emitting a photon. This is believed to be the first 
Simultaneous observation of these two profiles in the dayglow. 
b. Adjustments to Franck-Condon Factors 
The first set of spectra to be fitted to the data included 


Nitric Oxide y, 6, and € bands, Nitrogen LBH and V-K bands, OI 
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Figure 5-13 Intensity plots for LBH & VK showing quenching 
of VK. 
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2972A, OII 2470A, and OIII 2853A emissions. The first fine- 
tuning of the fit was accomplished by plotting the fit and 
data simultaneously in 400A segments and looking at each band 
for trends in undershooting or overshooting at a particular 
altitude, or for all altitudes at a particular wavelength. 
Any trend was then adjusted up or down in the original model 
spectrum, a procedure equivalent to varying the Franck-Condon 
factor (the factor which quantifies probability of transition 
for the synthetic models). This process was very tedious and 
time intensive, but a second fitting yielded an average 
decrease of 30% in y%. 

After the first set of changes had been made to the 
Franck-Condon factors, the same process of plotting and 
observing trends was repeated. Several unidentified peaks 
were seen in the data which are believed to be unidentified 
atomic emissions, Since their profiles matched that of an 
impulse function convolved with the instrument slit function. 
Other than the well-documented NII 2143A doublet, which was 
added as a synthetic spectrum, these emissions had not been 
previously reported in the literature, so the "Other Atomics" 
Spectrum was created (see discussion above). The two 
parameters, NII and "Other", were included in the fitting 
routine. Once these parameters were included, the adjustments 
to the Franck-Condon factors found for the Nitric Oxide bands’ 
synthetic spectra confirmed those used by Mack (1991) and 


Clayton (1990). The adjustments to the Franck-Condon factors 
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are summarized in Table 5-3. These changes to the accepted 
probabilities of transition are believed to be actual 
corrections and not just a function of the instrumentation. 
It 1s expected that further independent measurements will 
verify these corrections. 

The intensity profiles were derived from the final 
parameters of the fitting routine, but analysis of the data 
includes looking not only at the profiles but at the fits 
themselves. A complete set of plots of the data showing a 
comparison of data and fits, as well as accompanying plots 
showing the breakdown of the fits by constituent, are in 
Appendix E. Figures 5-14 and 5-15 show sample fits of the 


entire wavelength range at low altitude and high altitude. It 
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ails 


is apparent from these figures that the overall fit was 
excellent. 

It should be pointed out that there is a disagreement 
between the data and fit at the wavelengths where the nitrogen 
Second Positive emissions are present. Additionally, the only 
other poorly fit feature was the 2930A to 2990A range at 
altitudes above 200 km. This is one of the areas one would 
imeem eo bewable to fit well” since it includes the Of 2972A 
emission. Unfortunately, two of the dropped data points are 
at 2969.7A and 2972.9A, exactly where the peak of OI is 
expected. Also, there seems to be another emission at about 
2955A which is not explained by the VK bands or OI emission, 
creating three peaks in the data where theory predicts only 


two. This feature is shown in Figure 5-16. 
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Figure 5-16 Averaged spectrum for 285km. The peak at 2940A 
is VK(0,7), the peak at 2980A is due to OI 2972A, and the 
feature at 2955A is unidentified. 
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c. Temperature Profile 

In addition to the column density profile, intensity 
profiles, and fits of model spectra to data, a temperature 
profile also was determined. The temperatures described by 
Mack (1991) were used as initial guesses. Several of the 
temperatures at the lowest altitudes were increased during 
fitting to improve the match of synthetic spectra to data. 
The resultant temperature profile is shown in Figure 5-17. 
The temperatures determined by the fitting procedure are all 
considerably hotter than predicted by MSIS. In particular, 
the asymptotic temperature is 400° hotter. This is possibly 
due to a magnetic storm which occurred about three hours prior 


to launch. 
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Figure 5-17 The temperatures used by the analysis are 
indicated by symbols. The solid line is the MSIS 
temperature profile for the day of the launch. 
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VI. CONCLUSIONS 


A. SUMMARY OF FINDINGS 

Thermospheric dayglow data from a rocket-borne 
spectrograph were analyzed in this7~ study. Computer 
synthesized spectra were created for the N, Lyman-Birge- 
Hopfield and Vegard-Kaplan bands; NO y, 6, and € bands; OI 
2972A, OII 2470A, and OIII 2853A atomic lines; and NII 2143A 
doublet feature. A spectrum of other possible atomic 
emissions and a background factor were also considered in the 
fit. These synthetic models were fitted to the observed data, 
then the scaling factors were used to determine an intensity 
profile (column density profile in the case of nitric oxide). 

The intensity and column density profiles obtained were 
within 10-20% of the values acquired by simpler fitting 
routines in all cases except the OI 2972A emission. Thus the 
analysis supports the concept of fitting of just one feature 
to arrive at density and intensity profiles. A relationship 
was noted between VK and LBH emisSions that would seem to 
support the theory of quenching of the nitrogen VK system. 
The data above roughly 3000A could not be modelled without a 
N, Second Positive synthetic spectrum. The loss of datapoints 
was a major detriment in fitting the OI 2972A emission, as 


there were two data dropouts where a peak was expected. The 


aye) 


temperature profile which best fit the data was higher at low 


altitudes than previous analysis had predicted. 


B. TOPICS FOR FURTHER INVESTIGATION 

With another launch scheduled for March 1992, further 
Spectral analysis of the March 1990 data is imprudent. 
Contemplating the prospective data, a synthetic spectrum for 
the N, Second Positive would enhance the analysis. It is 
expected that there will be no problem with data dropouts on 
the next flight (Quint, 1991), so a repeat of this analysis 
(with the new data) will greatly improve the accuracy of 
fitting of the OI emission and possibly show that it can be 
accurately fitted over a limited wavelength range. 

The quenching of VK below 200km should be paralleled by 
the density of the quenching agent, aS well as a rise in 
emissions corresponding to the relaxation of the quenching 
agent from its excited state. The quenching phenomenon showed 
promise for rewarding study, but was not the subject of this 
thesis. 

Looking toward the future, the MUSTANG will be part of an 
operational satellite in the late 1990’s, providing a constant 
flow of data for analysis. Considering the numerous new 
mathematical theories for computer network calculations, there 
may be a method better suited for fitting large quantities of 
data than the GRIDLS routine. A modification to the fitting 


programs to place less stringent requirements on the initial 
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programs to place less stringent requirements on the initial 
guess and to provide a more formal check for global minimum 
would greatly improve the effectiveness of the existing 
routine. Additionally, automation of the analysis, including 
derivation of the electron density profile, would certainly 
add to the viability of MUSTANG as a provider of geophysical 


parameters for national defense agencies. 
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APPENDIX A 


The following table shows the unshifted and shifted array 
elements of data, and their corresponding wavelengths. 
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APPENDIX B 


Appendix B contains the least squares fitting program and 
its subroutines. Additionally, the synthetic spectrum 


generation routines are given. 


Walden, 
@getfunc 
@fchisqr 
pro nimp 


a 


“=e “Se “Se “=e Se Ne Se We Be Ve We Be We We We We Se BWe We We Wo We We We We We 


Too, 


fit oie 


modified from Mack,1991 


This program fits data for the 1800 to 3400A range with models of the 


NO gamma, 


delta, 
oxygen emissions at 24/70, 


Kaplan bands, 


at 2143A, 
through an iterative process. 
fchisgr.pro as subprograms; 
and it writes the fit by altitude to nimpfit all.dat. 


various atomic emissions, 


and epsilon bands, N2 LBH emisSions, N2 Vegard- 

2972, & 2853A, N+ doublet 

and the background level 

Tt uses.gridls.pro,.getfunc. pro ;eana 

it useS parameter.dat as an in/out file; 
The subprograms 


utilize synthetic spectra which are stored in allvk.dat, all _lbh.dat, 
nii 2143doub.dat, nodlt.dat, nogammack.dat, noeps.dat, oiii 2853.dat 
olin 2470.dat, Oae29/2. dat, andwether spec jac 


parameter.dat - file of parameters- records have 25 elements each 
record#— contents 
O - temperature 
1 - altitude record number (1 through 25;1-17 up,18-25 down) 
2 - chi squared of the fit 
3 - scale factor for O 2972A emission, in Rayleighs(intensity) 
4 - sigma for O scale factor 
5 - scale factor for O+ 2470A emission, in Rayleighs 
6 - sigma for O+ scale factor 
7 —- scale factor for O++ 2853A emission, in Rayleighs 
8 - sigma for O++ scale factor 
9 - scale factor for the N2 LBH bands 
10 - sigma for LBH scale factor 
11 - scale factor for the NO gamma, delta, and epsilon bands, in 
molecules per square centimeter (column density) 
12 - sigma for NO scale factor 
13 - scale factor for the N2 V-K bands 
14 - sigma for VK scale factor 
15 - scale factor for N+ 2143A emission, in Rayleighs 
16 - sigma for N+ doublet scale factor 
17 - scale factor for synthetic atomics spectrum, in Rayleighs 
18 - sigma for atomic spectrum 
19 - scale factor for the background, in Rayleighs per Angstrom 
20 - sigma for background scale factor 
21 - altitudes in kilometers 
para - array((2*nterms)+2) in same order as record order in file 
a - array(nterms+l1) of temp followed by scale factors 
yfit - 512 element array of best fit 
start- first record number 
stop -— last record number 


openr,30,’(mustang]m90 int.dat’ 


;open data file 


aa=assoc(30,fltarr(512)) 


openu, 31,  nimptiitgaliee dae. 


;open file to save yfit 


bb=assoc(31,fltarr(512)) 


openu,32,’parameter.dat’ 


;open parameter file 


PPp=assoct 326 lta mite oe 


openr, 33,’(walden.datafix]m90 domwt.dat’ 
dd=assoc(33,£fltarr(512)} 
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Sleshe eee) 


stop=511 
npts=(stop-start)+1 Pi Oe POINntswwompe fit 
nterms = 9 ;# of terms to consider, change if 


;parameters are added 
para=fltarr(2*nterms+3) 
a=fltarr(nterms+l ) 


for i=0,24 do begin 


for j=0,(2*nterms)+2 do begin ;£i111 para with values for alt i 
hold=pp(j) 
para(j)=hold(1) 

endfor 

alt=para(1) 

count=0 

a(0)=para(0) ;assign first element in a temp 

for j=3,(2*nterms)+1,+2 do begin “fiil aiweth scale factors 


count=count+l 

a(count)=para(j) 
endfor 
y=aa(start:stop,alt) 


; 
;deltaa is step size for estimations, cannot be zero, first element is dummy 


h=a 
for j=l,nterms do begin 
h(j)=a(j) 
if (h(j3) eq 0.0) then h(3)=1000.00 
endfor 


deltaa=h/10. 


mode=1 ;use sdom for weighting 
weights=dd(alt) 


gridls,y,weights,npts,nterms,mode,a,deltaa,sigmaa,yfit,chisgqr,start,stop 
bb(alt)=yfit 


count=0 ;Ssave parameters from best fit 
para(2)=chisgqr 
for j=3,(2*nterms)+1,+2 do begin 
count=count+l 
para(j)=a(count) 
para(j+1)=sigmaa(count) 
endfor 
for 3=0,(2*nterms)+2 do begin ;put parameters in appropriate 
hold=pp(j3) ;record in parameter.dat 
hold(i)=para(j) 
pp(j)=hold 
endfor 


Bent, alt=',alt,’a=',a,'’chisqr=’,chisqr 


| endfor 
close,30 
close,3l 

| eLose, 32 
close, 33 

return 

end 

@gridls 
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*Walden,1991 modified from Mack,1991 
Puncrion Getrume, a 


: This function computes the synthetic spectrum from 1800 to 3400 A 
; for given scaling factors. 
: INPUTS 

; a - 9-D array of parameters 

; a(Q0) - temperature at altitude of interest 

; a(l) - scale factor for O emission @ 2972 
: a(2) - scale factor for O+ emission @ 2470 
> a(3) - scale factor for O++ emission @ 2853 
' a(4) - scale factor for N2 LBH bands 
; a(5) - scale factor for the sum of NO gamma, delta, 
; and epsilon bands 
; a(6) - scale factor for N2 Vegard-Kaplan bands 
; a(7) - scale factor for N2 doublet emission @2143 
; a(8) - scale factor for "Other Atomics" spectrum 
; a(9) - scale factor for background : 


GETfunc returns an array of 511 points to be compared to the 
experimentally recorded profile from 1800 to 3400 A. 


LBH - The N2 LBH bands convolved with the slit function 
NOSYN - The sum of the NO gamma, delta, and epsilon bands convolved 
with the slit function 


=s “~e sc We “oe Be Ne We We We We We 


VK - The N2 V-K bands convolved with the slit function 

OI - The OI emission @ 2972 A convolved with the slit function 
Ome - The OII emission @ 2470 A convolved with the slit function 
OIII - The OIII emission @ 2853 A convolved with the slit function 
N2 - The N2 emission €2143 A convolved with the slit function 
AT - The spectrum of possible atomic emissions 


openr, 207 Bibeorieb a 
rr=assoc(Z0,7rttarei sl) 
LBH=rr(a(0)/100) 


openr,21, ’NOGAMMACK. DAT’ 
Ss=assoc(Zl ,fltamei 512) 
NOGAM=ss(a(0)/100) 


openr ,22, NObLL DAL 
tt=assoc( ZZ; 8learerol)) 
NODLT=tt(a(G)7 Fee, 


openr,23,’NOEPS.DAT’ 
uu=assoc(Z3, ft lLtare(siZ)) 
NOEPS=uu(a(0)/100) 


openr,24,’ALLVK. DAT’ 
vv=assoc(24,fltarr(512) ) 
VK=vv(a(0)/100) 


open, 25, Olmg72. Dae 
ww=assoc(25,fltarr(512)) 
OI=ww(0) 


openr, 26; “Ol 2470. Pal 
xxwassoc(26,fltarr(512)) 
Ora (ae) 


openr, 27, Ollie cosee ae 
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yy=assoc(27,fltarr(512) ) 
OIII=yy(0) 


openr,28,’NII_ 2143DOUB.DAT’ 
g@2=assoc(28,fltarr(512) ) 
N2=22(1) 


Spent,29,'otherspec.dat’ 
mesp=assoc(29,fltarr(512)) 
AT=atsp(1) 


NOSYN = NOGAM + .25*NODLT + NOEPS ;25% flourescent efficiency of 
s;delta band 


getfunc = a(4)*LBH + a(5)*NOSYN/3.133e6 + a(6)*VK + a(1)*OI §$ 
+ a(2)*OII + a(3)*OIII + a(7)*N2 +a(8)*AT + a(9)/3.133 
;background is divided by wavelength bin so that 
;scale factor is absolute intensity, NOSYN is divided 
;by 3.133e6 for same reason 


close,20 
close,21 
close,22 
close,23 
close,24 
close,25 
close,26 
Grose, 2/ 
close,28 
close,29 


return,getfunc 
end 


ay 


;Mack, 199" 
function fchis@nr,y,welghts ,npts, tree ;modenvi mens tare Uo 
This function evaluates the reduced chi squared 
for a fit te themaaecar 
Adapted from Bevington, pg 194. 
INPUTS 
Y -~ array of data points 
weights- array of weightings based on sdom’s of data points 


npts - number of data points being fit 
free - number of degrees of freedom 
mode —- determines method of weighting the fit 


+l instrumental weighting(sdom weights) 
O no weighting 
-1 statistical weighting 


™e “we We Se We We BWe We We We We We We We We 


VEVE - array of fitted points 

start - first record element to be fit 
stop - last record element to be fit 
count=0 

chisqr=0. 


if mode eq 0 then wt=1. 
if mode gt 0 then wt=weights 
if mode lt 0 then wt=1./abs(y) 


for i=start,stop do begin 
chisqr=chisqr + wt(i)*(y(count)-yfit(count))°2 
count=count + 1 

endfor 

fchisqr=chisqr/free 

bOCULN, © earache. 
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*Walden,1991 
aeeeeeoncalc,wl,int,slit ;Calculates synthetic lbh spectra at temps 
Spenu,10,’all lbh.dat’ ; From 100 to 2000 deg for wavelengths from 
jmem—assoc(l0,fltarr(512)) ; 1800 to 3400 A. Results are convolved with 
; MUSTANG instrument slit function. 


wl=indgen(512)*3.133 + 1798 


openr,11,’[walden]slit.dat’ 
Seeassoc(l],fltarr(l1l)) 





Set=s] (1) 
ie) 
fomee—!O0,2000,100 do begin 
i1=i+1 ; Wavelength 
int=wl*0 ; variables for tweaking individual transitions 
intl=wl*0 
int2=wl*0 
int3=wl*0 
int4=wl*0 
ima ow 1 * 0 
int6=wl*0 
int7=wl*0 
int8=wl*0 
int9=wl *0 
intl0=wl*0 
lbhsyn,t,4,9,wl,int hs @ eG) 
Mensyn,t,9,6,wl,int peace 4 .6 
memeyn,t,5,10,wl,int 7 lst 7 9 
menmsyn,t,1,7,wl,intlod pO e:.-6 
int=int+intlO*1.7 
Monsyn,t,6,ll,wl,int oS 3a 0 
Mensyn,t,2,8,wl,int eee 22 
Pensyn,t,3,9,wl,intl 7 155378 
int=int+intl*1.203 
memsyn,t,4,10,wl,int ceo OG 
Mensvyn,t,0,/7,wl,int hs) Ieee 
iensyn,t,5,l1,wl,int2 > 1887.9 
Pit —int+int2*1.4 
teasyn,t,1,8,wl,int3 elogea. 2 
Int=int+int3*1.3 
Measyn,t,6,12,wl,int 7 90S ..4 
memsyn,t,2,9,wl,int > 1910.9 
mensyn,t,3,10,wl,int4 eee 1S 
int=int+int4*.8 
mensyn,t,4,11,wl,ints elo 450) 
int=int+int5*.9 
Mensyn,t,0,8,wl,int eee 1 S)Sy Ses, 
Mensyn,t,5,12,wl,int6 » Prge2.4 
int=int+int6*.757 
| memsyvn,t,1,9,wl,int pe lO 2a 
| imsvn,t,6,13,wl,int O80. 1 
Monsyn,t,2,10,wl,int ; 1989.6 
| Mensyn,t,3,1l1,wl,int? 2 OO eG 
int=int+int7*.8 
| ionsyn,t,4,12,wl,ints ete 2a 2 
int=int+int8*.9 
ieiicvin,t,0,9,wl,int OAS 
ioisvn,t,5,l3,wl,int9 > 2041.9 
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sgher nonce se gye sli 


lbhsyn,t3 1,007 wie ane > #820568 
lbhsyn,t,6,14 wl, ine 23 2059. 
lbnsyn,t,2, 2 ew le = 209 3: 
lbhsyn,t,; 37 2a wea 2 20 12 
LDASYN, 747 b37v le 2 US = 
lbhsyn,t, 5) 4a. ot PAILS) 
lbhsyn,t,0, 1077) ine ye Ne 0, 
Ibhsyn, t ,6 25), wlan 2 a ae 
lbhsyn , ee wale ines oe 2b da, 
Ibhsyn,t,2, 12, wee 208 SG 
lbhsyn,t,3,13, wae * we 2 ees 
lbhsyn,t,4,14,wl,int Re PAA 35 
lbhsyn,€,5,157 70 ue ae 2 ee, 
LeASyYN;t 70k wee Hee oh rT 
bhasyn,t,6, 16, wine 2 eZ235 
Pohs yn, lk, 125 ween mn 2a 
SiS yiyee, 2, L3;,WL, ine patie PA(enl 
ibhisyn, ts 3,04, wl pane ee A, FS 
Tbhasyn,t, 4, 153.070, 1nG - 2297 
EDS Vint, oy loi lem Zoe 
Pbhsyn, & 6107 jw one Pez 3 a2 
TBHS Yn; Ey Lp See eit ; 2349 
lbhsynjt, 2,04 wil ae 2 Eels 
IL DAS YNe C73 7S Lee - 2304 
lbhsyn, 747 lo yw lee a 20 
lbhsyn, toys ee °* 2419 
lbhsyn, € Oyo wlan c > 92437 
Iibhsyn, € 72,457 w ieee PAA TES 
IDNs VN tps, lew, pe >; 2496 
Ibhsyn,b,471 wl ane Zona 
IpbASYVn7 t,o, fo, wl ae oe Seoul 
lbhsyn,t;6,19; 7) poe s7- 2549 
lbhsyn,t;4, 06 .wl) mt me AS SIG) 
IbhAsyn,€, 2,278 ee AL be 25 oS. 
lbhsyn,t;6,20; wil. int 222690 
lbhsyn,t;76,21,"), ine - 2801 


Print; total~inr, 


spec=convol(int,slit) 
lbh(i)=spec/total (spec) 


endfor 


close,10 
elosevim 


eturn 
nd 


[SYNSPEC]LBHSYN 


© ~~] CON © OO?’ 


~l Fe 
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swalden,1991 
feemalivkcalc,wl,int,slit 
openr,1,’[WALDEN]slit.dat’ 
ag@enu,s,’allvk.dat’ 
sl=assoc(1,fltarr(ii)) 
slit=sl(i) 
Bemetescoc(3,fltarr(512)) 
vl=indgen(512)*3.133 +1798 


L= 


For t=100,2000,100 do begin 


1=i4¢+l 


int=wl*0 

intl=wl*0 
int2=wl*0 
int3=wl*0 
int4=wl*0 
int5=wl*0 
int6=wl*0 
mc 7=wl*0 
int8=wl*0 
int9=wl*0 

ite tL O=wl*0 
intll=wl*0 
intl2=wl*0 
watd 3=wl*0 
intl4=wl*0 
ome) 5=wl*0 
intl6=wl*0 
Mae) /=wl*0 
int18=wl*0 
intl9=wl*0 
mae 2 0=wl *0 
mates l=wl *0 
int22=wl*0 
iat 2 3=wl*0 
mat 24=wl *0 
mee >5=wl *0 
ime zo=wl*0 
int27=wl*0 
int28=wl*0 


This section calculates the transitions’ band systems, summing each one 


Meco the total spectrum. 


TRANSITION 


meoywn,t,6,1,wl,int 
mesyn,t,4,0,wl,int 
een ,t,13,5,wl,int 
meosyn,t,11,4,wl,int 
eesyi,t,9,3,wl,int 
mesvil,t,7,2,wl,int 
mesvn,t,5,l,wl,int 
Msyn,t,3,0,wl,int 
miesyn,t,12,5,wl,int 
mmsvn,t,10,4,wl,int 
wesyn,t,8,3,wl,int 
maesyn,t,6,2,wl,int 
wieewt,t,4,l,wl,int 
Pinoy, t.3,60,Wl,int 


;Calculates synthetic Vegard-Kaplan 

mee PecCtra -aumcemps Erom 100 to 2000 

; degrees for wavelengths from 1800 to 

; 3400 A. Results convolved with MUSTANG 
; instrument slit function 


; these variables provide a means of tweaking 
; Imaavidual transitions 


WL 


ALIS 
1808 
PS23 
e224 
S27 
1833 
1841 
siebeye: 
1860 
ieilere: 
1868 
1876 
A ehis\ 7 
see S) 
8] 


™me “se ®e ~se We WO WE We We We We WE We WE 


VKESVny t,) bp ine 
vksyn, &, 2, 0, wu 
VKSYVIl, Gao 7, dye 
VESVMyet sot oy en 
VeSYN, Cy 5727 WL me 
VkKsyn,t, 3, 4,0) me 
vVksyn,t, 2,657.9, 
VeSyvn;,t, 10757 Ween 
VESVI,© > Oy 47, Ww le mieniite 
vVksyn,t,1,0, ie 
IME=INC+ Ine Ls Peo 
VKSV¥N; tC, 6737p 
VKSYN ,-t7 4) 2, wie 
VKSVn;7t,13 57, wine 
vVesvn, t; ll, 6,v7 ine 
vksyn,t,9;o,wl,2nt 
WS Vil, © ye gly we 
VASVN, 6,17, 4, We 
Veeyil, ey 070, WL elm 
int=int+int2*2.6 
WieSVTl 7:6 po posal 2, mente 
VKSVhyty lle pw pane 
VKESVN; t;3;2 Vl, Ine 
vViSvn, t,l0,6,W ime 
VESVN, 6, Oponw ly Lime 
VESYN, t > Jee, wees 
int=int+int3*1.5 
Vkoyn,€, 67470 pie 
VESVI; 0,475 Wie me 
VESVn, ©) loc, viene 
VESVly oy Lie weet 
VK SVaiG eG, 90,6 ete 
VKSVNGt, 2, 2,Wl, on 
vkSyn,t, 7) pow eae 
VKSYN, €, S94 Ww lee 
VRSivIy ec 0 nw eee 
int=int+int4*2.6 
vVksyn,t;lZ76)W ame 
VKSVIN-t, LUe | wil rit 
VESVI) Gy oy SW lene 
VESVN » eyo, 6, Wil ane 
VESVN, Cp Oy oy Wl pene 
VieSV DiseeieZ | Wile. 
int=aincrinto® lo 
vVksyn,t,13,9,wl; int 
vksyn,t,4,4,wl,int 
VKESViIne ct, Lis iim 
VKSyny t, 959750, ne 
Vi ovnit, 25 poe 
VES VIN, Ce, ) Op Wwl nt 
VESVIi pty oy Ole ie 
vksyn,t, /2,0, 4) ant 
VESVRy G02, eee 
int=int+int6*2.6 
vksyn,t, 1076, 71, ane 
vksyn,t,3,4,wl,ane 
VKSYN to op ay, Ie 
vksyn,t, 6, 6,W.L ine 
VESYN pty saw fe 
int=int+tint)- 25> 


=e “se Se Seo “6 Se We we We =e “oe “e Se De We =e “=e “oe Se Se Se Se We =e “se Se Se Se We Be We Se We 


™=e “se “~e “We We Se We We We "=e “*®we Se “Se Se WE 


=e “Se “Se Be WO 


1900 
1901 
ILS) 05) 
iy 
1923 
Pos 6 
Ihe) 253) 
1943 
1350 
1954 


Ine e 
eS 
1960 
196 3 
1969 
IS 
ISSN 
2010 


PALO }A tsa 
2024 
ZZ 
Z02Z9 
2038 
2047 


2050 
Z0G> 
2067 
Z07 1 
2079 
Z0iS Ss 
2090 
Z06 
2109 


BIAS 
2122 
22> 
Zi3s 
2147 
2148 


753 (5y 
2166 
2167 
Zale 
2189 
2191 
2209 
2213 
PAP ad 


Z2ce 
22351 
2236 
ZS 5 
2250 
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mesyvn,t,13,10,wl,int 
Mesyn,t,11,9,wl,int 
mrosyn,t,4,5,wl,int 
mesyn,t,9,8,wl,int 
Prescyn,t,7,/7,wl,int 
wesyn,t,2,4,wl,int 
feeeyn,t,12,10,wl,int 
mesyn,t,5,6,wl,int 
mroyn,t,10,9,wl,int 
fPeooyin,t,0,3,wl,int8 
mime=—int+int8*2.6 
mesvn,t,3,5,wl,int 
mesvyn,t,8,8,wl,int 
mesvn,t,6,7,wl,int 
masyn,t,13,l1,wl,int 
meoyn,t,l,4,wl,int9 
fme=int+int9*1.5 
waeoyn,t,l1,10,wl,int 
vksyn,t,4,6,wl,int 
mesvyn,t,9,9,wl,int 
vksyn,t,7,8,wl,int 
meeyn,t,2,5,wl,int 
asym, t,l12,11,wl,int 
mesyn,t,5,/7,wl,int 
feeoynm,t,10,10,wl,int 
vksyn,t,0,4,wl,intl10 
mmt—int+intlO*«2.6 
meoyn,t,8,9,wl,int 
Mreoyn,t,3,6,wl,int 
feeeyn,t,13,12,wl,int 
feeeyn,t,6,8,wl,int 
fee yn,t,l1,11,wl,int 
fees, t,l1,5,wl,intlil 
mie—1ntt+intli*1.5 
Mimoyn,t,9,10,wl,int 
mesyn,t,4,/,wl,int 
Mmoyn,t,/,9,wl,int 
meovn,t,2,6,wl,int 
tmecyn,t,12,12,wl,int 
mesyn,t,5,8,wl,int 
meeyvn,t,lO,ll1,wl,int 
meoyn,t,8,10,wl,int 
fiesyn,t,0,5,wl,intl2 
mit —=int+intl2*2 .6 
mesoyn,t,3,/7,wl,int 
fioyn,t,13,13,wl,int 
mroyn,t,6,9,wl,intl3 
mie—iLnt+intl3*.83 
weovn,t,l1,12,wl,int 
meoyn,t,L,6,wl,int14 
Pme=1nt+intl4*1.5 
Mmaoyn,t,9,ll1,wl,int 
wesyn,t,4,8,wl,intls 
mae=iInt+intls*.67 


feesoym,t,7,l0,wl,intl6 


int=int+int1l6*.6 
mecyn,t&,12,13,wl,int 
weSsym,t,2,/,wl,int 
meoyn,e,l0,12,wl,int 
evn ,t75,9,wl,int 


™=e we =e BZe BW WE We WE WE =e =e Be BE WE ™=e “se =@ Be Be We We We WE WS 


™m™e se =6¢ Be WE WE 


™me Me Oe We We WE We We WE 


me We We WE 


2202 
220 
22S 
PSII 
2209 
2302 
Dn SP, 
Zea 
Zio 
2334 


2347 
2347 
2306 
2302 
239 


2382 
2394 
2397 
2417 
2425 
2435 
2443 
2449 
2463 


2469 
2474 
2490 
2494 
2504 
Zoe 


2522 
2524 
2547) 
2508 
Zon 
Zotd 
2576 
2602 
2605 


ZG a3 
7a\6)7al\) 
Zo 32 


ZOS 7 
Zoe 


2659 
2668 


2689 


ZOg56 

Ze 

ISS, 

2724 
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VKSVMy t , Oe len 
vVksyn,t,13, f¢7 0 int 
vksyn,t, 0,6, wily lL? 
int=int+intl7*2.6 
vksyn,t,3,8,wl,int 
vksyn,t,117 13 wines 
int=int+intl8*.7 
vksyn,£,6,.0, wl ine 
VKSV ily by 9 pew ee 
Vksyn;t; ly 17 Wile 
int=int+intl19*1.5 
vVksyn; t,4;, 9 .w lnc 
visyn, ct, 7, ll wa mt 
vksyn,t,12,1¢,wi,ineZo 
int=int+int20*.67 
visyn,t, 10,03 wi, ine 
Vikevn tic, 6, Wioermes | 
INnt=inttint2Zl= 67 
Visvim, t, > 710 wiles eZ 2 
int=int+int22*.5 
vksyn,t, 0,102, wig 1nt23 
int=int+int23*.5 
VINSYN, © 13, Loewen inet 
VKSVYN,t,0,/ pw mes 4 
int=int+int24*2.6 
vksyn,t,3,9,wl,int 
VKSYN pte wl ne 
vksyn,t,6,11,Wwieaine 
VKESYN E> 7, Loe ae 
vVisyn,t, 1,8; wi tatZ oO 
int=int+int25*1.5 
VKSyVn,t;,4,)0, Woe mt 
VkSsyn;t;l2,15 want 
VKSYN, tC, ipl ZV leant 
VKSYN,t,10, 14 wi int 
VKESYN ty 2 wile, ome 
VESVN, C7) opt wie Tit 
viegyn, tls Gy vee nt 
VRSYN;t so, 13 went 
VKSyn) tell, lSmwiein 
vksyn,t,3,10,whant 
VEsyn,t;,07e," ll, amt 26 
int=int+int26*2.6 
VKSYN,t;0, 1.2.0 pene 
Vksyn,t;9, 14,02, int 
VKsyn,t, 12,167,027, nc 
vksyn,t,4,11)wiaint 
vksyn,t,;1, 9, whine 7 
Int=1nNt+ines / “a> 
VKSVN> 6597 7s) wien & 
vksyn,€,107 ow ent 
VKSVMjt, 27 0 in © 
VESYN, Ee, Ona ween 
vksyn, t, le, awl ane 
VK SVN pty 6 lo Wile © 
vksyn,t,;11,16 ; whe int 
VKESVM, ty 31 1 ween © 
VESVNjt, Oop lee 
vksyn,t,0,9,.2,int2Zs 
int=int+int28*2.6 


™~e S~e =e Be WE =e %e %e Re Ze Be We We We BW WE =e %»ec %e Re WE 


=e *®e Be Be Be We Be We BWe We 


2748 
2761 
Za 62 


Za 7 
2762 


2783 
2810 
2818 


2826 
2845 
2848 


ZB1S 
Zo 1d 


Aes | 
2909 


Zo 6 
2937 


Z3'8 
2942 
74 {8) 3 Ah 
I, NIN 
Zoo G 


3002 
3013 
S017 
3046 
3062 
3068 
3087 
3087 
B19 
3128 
S33 


S43 7] 
e1oSis, 
SIEM 
g29 8 
e200 


S20 
3235 
2) PASH, 
5270 
3276 
3285 
5) ons: 
3342 
3345 
3353 
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miesyl,t,9,15,wl,int 
Pesyn,t,12,17,wl,int 


wari) -convol(int,slit) 


emaror 

close,l 

close, 3 

return 

end 
@[synspec]vksyn 
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»;Walden,1991 
PRO NODLTGALE ,WLE,INT,cobrLl 


OPENU,10,’NODLT.DAT’ ;This procedure calculates synthetic 
DLT=ASSO0C(10,FLTARR( 512) ) ; Nitric oxide delta band spectra 

; at temps from 100 to 2000 deg and 

; convolves results with MUSTANG instru- 

- Ment slit function. 
WL=INDGEN(512)*3.133 + 1798 


OPENR? 11 7’ (WALDENISEIT DAT’ 
SL=ASSOC(11,FLTARR(11) ) 


SLIT=SL(1) 

I=0 

FOR T=100,2000,100 DO BEGIN 
I=I+1 ; Wavelength 
INT=WL*0 
DEBSYN 7 T On G, Why oNT ; 1909.8 
DEES 7 20 aioe ry dS tSh OURS 
DEGSYN, F,.0;72 bk, INT ames liv aera’, 
DEESYN, 2, 0 p37 tN 7.2854 20 
DELSYN, T,0,;4,Wh,i1NT a 2 Sa 
DEGSYN, T7045), Wie ee * i235 Ol 
DEESYN, T3067 4 by en 2 . 2409 46 
DELSYN, 1,07 77 WE UNE ; 2509 
DELSYN, T7005 8307 1 ths 412 P48) 
DELSYN, 2,0) 93 Vibeeher , 92739 
DELSYN, TU 00 Wea T s Zoo 
DELSYN, TUR WeanwlNE + he 003 
DELSYN, TU; 12, WG, EN? ot SpE OZ 
DEESYN, 1, 0; U3 Wb, nr prmeo3 
DEESYN; 17 07 4 Wien ; 3488 


DLT(I)=CONVOL(INT,SLIT) 
ENDFOR 


CLOSE VLG 
CLOSE 


2ETURN 
IND 
i[SYNSPEC]DELSYN 
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Walden,1991 

This program derives the synthetic spectrum using tweaking factors 

first determined by Mack(1991) for 1800-2100A, and extended the full range 
of the MUSTANG data. 

RO NOGAMCALC,WL,INT,SLIT 


OPENU,10,’NOGAM.DAT’ ;This procedure calculates synthetic 
GAM=ASS0C(10,FLTARR(512) ) ; nitric oxide gamma band spectra 
Pe-TNDGEN(512)*3.133 + 1798 ; at temps from 100 to 2000 deg and 

; convolves results with MUSTANG instru- 
> ment slit function. 


Seer, tl,’ {(WALDEN]SLIT.DAT’ 
SL=ASSOC(11,FLTARR(11) ) 
SLIT=SL(1) 


I=0 
FOR T=100,2000,100 DO BEGIN 

I=I+1] 
;These variables allow for the adjusting of 
Menagaividual transition’s probability 
INT=WL*0 
INT1=WL*0 
INT2Z=WL*0 
INT3=WL*0 
INT4=WL*0 
INT5=WL* 0 
INT6=WL*0 
INT/=WL*0 
INT8=WL*0 
INT9=WL*0 
INTLO=WL*0. 
INTL1I=WL*0. 
INT12=WL*0. 
INT13=WL*0. 
INT14=WL*0. 
INT15=WL*0. 
INT16=WL*0. 
INT1L7=WL*0. 
INT1L8=WL*0. 
INTLO=WL*0. 
INT2Z0=WL*0. 

| INT2Z1=WL*0. 

) INT22=WL*0. 

| INT23=WL*0. 

| This section calculates the spectrum from each transition’s band 

system. 


Transition WL 
GAMSYN,T,3,0,WL,INT ee 9 56; 1) 
GAMSYN,T,3,1,WL,INT meZ029. 6 
GAMSYN,T,2,0,WL,INT9 7 2046..5 
INT=INT+INT9*.75 
GAaMsYN,T,3,2,WL,INT ene > 
GAMSYN,T,2,1,WL,INT ae? leo 4 
GAMSYN,T,1,0,WL,INT1 ree Aree! 
PNT=INT+INTL*1.125 
GAMSYN,T,3,3,WL,INT > 2192.8 
GAMSYN,T,2,2,WL,INT2 > 2215.4 
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INT=INT+INT2*1.26 
GAMSYN,T,1,1,WL, INT3 
INT=INT+INT3*1.3 
GAMSYN,T,0,0,WL,INT4 
INT=INT+INT4*1.2 
GAMSYN,T,3,4,WL,INT 
GAMSYN,T,2,3,WL,INTS 
INT=INT+INT5*1.6 
GAMSYN,T,1,2,"l, ENGI 
INT=INT+INT10*1.8 
GAMSYN,T,0,1,WL, INT6 
INT=INT+INT6*1.2 
GAMSYN,T,3,5,WL,INT 
GAMSYN,T,2,4,WL, INT 
GAMSYN,T,1,3,WL,INT7 
INT=INT+INT7*1.6 
GAMSYN,T,0,2,WL,INT8 
INT=INT+INT8*1.2 
GAMSYN,T,3,6,WL,INT 
GAMSYN,T,2,5,WL,INT 
GAMSYN,T,1,4,WL,INT17 
INT=INT+INT17*1.3 
GAMSYN,T,0,3,WL,INT11 
INT=INT+INT11*1.2 
GAMSYN,T,3,7,WL,INT 
GAMSYN,T,2,6,WL,INT 
GAMSYN,T,1,5,WL,INT12 
INT=INT+INT12*1.3 
GAMSYN,T,3,8,WL,INT 
GAMSYN,T,0,4,WL,INT13 
INT=INT+INT13*1.2 
GAMSYN,T,2,7,WL,INT 
GAMSYN,T,1,6,WL,INT14 
INT=INT+INT14*1.3 
GAMSYN,T,3,9,WL,INT 
GAMSYN,T,0,5,WL,INT15 
INT=INT+INT15*1.2 
GAMSYN,T,2,8,WL,INT 
GAMSYN,T,1,7,WL,INT16 
INT=INT+INT16*1.3 
GAMSYN,T,3,10,WL,INT 
GAMSYN,T,0,6,WL,INT18 
INT=INT+INT18*1.2 
GAMSYN,T,2,9,WL,INT 
GAMSYN,T,1,8,WL,INT19 
INT=INT+INT19*1.3 
GAMSYN,T,3,11,WL,INT 
GAMSYN,T,0,7,WL,INT20 
INT=INT+INT20*1.2 
GAMSYN,T,2,10,WL,INT 
GAMSYN,T,1,9,WL,INT21 
INT=INT+INT21*1.3 
GAMSYN,T,3,12,WL,INT 
GAMSYN,T,2,11,WL,INT 
GAMSYN,T,0,8,WL, INT22 
INT=INT+INT22*1.2 
GAMSYN,T,3,13,WL,INT 


GAMSYN,1T, 1720, Whe iniZs 


INT=SINTTINGZS4* 13 
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GAM(I)=CONVOL(INT,SLIT) 
ENDFOR 


CLOSE,10 
SboaE,11 


RETURN 
END 
@[ SYNSPEC ]GAMSYN 
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;Walden,1991 
PRO NOEPSCALC 


OPENU,;1, NOBES .PAT- 
EPS=ASSOC(1,FLTARR(512)) 


PIN procedure calculates synthetic 
nitric oxide epsilon band spectra 

at temps from 100 to 2000 deg and 
convolves results with MUSTANG instr- 
ment slit function, where EPS(1) 
corresponds to T=100,etc. It uses 
scale factors determined by Bosserman 
except for (0,1), which Bosserman 
used .84 to tweak 


™~e “*®*s “Xe Be De Ws We WE 


WL=INDGEN(512)*3.133 + 1798 


OPENR,2,’[([WALDEN]SLIT. DAT’ 
SL=ASSO0C(2,FLTARR(11) ) 


SLIT=SL(1) 

I=0 

FOR T=100,2000,100 DO BEGIN 
I=I+1 ; Wavelength 
INT=WL*0 
INTL=WL*0 
INT2=WL*0 
INT3=WL*0 
INT4=WL*0 
BP Sow Ne flO Wily. lie 7 ly oo 
EPSOYN? fez, 2) Wo, Le - 21349). 5 
EPSSYN? a, 1, Woe lie el silos 7 
ERS SMNe § 0e.0 . Wig laa oardo 1076 
INT=INT+INT1*.63 
EPSSYN; 2.2, 3,WL, om 2 oS 
EPSSYN il, > Wi DN pS a. 
EPSSYN, 7,0, 1,WL,2na2 > bo 45. 0 
INT=INT+INT2*.90 
EPSSYN{ fl, 2 VoL ; 1998.6 
EPSSYN, 2,072, 0b) Uae Fc OL as 
INT=INT+INT3*1.03 
BPSSYN; 12,4, Wh poe pee 30 
EPSSYN, £403 ,WL, nud > 2094.5 
INT=INT+INT4*1.11 
EbSovl ol, ie - S254 5 
EPSSYN? 10,4, WL, UNS aly Gee 
EBPoSYN, T7290, WL, Nm 7 t2294.0 
EPSSYN,T,2,9,WL,INT Pp 230550 


EPS(I)=CONVOL(INT,SLIT) 


ENDFOR 
CLOSE, 2 
CLOSE, L 
ETURN 
ND 


[synspec]EPSSYN 
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;Walden,1991 
moowoin27/2calc ;This program creates a synthetic model for a 
; neutral atomic oxygen emission 


Meenu,t,’OL 2972.dat’ 
dat=assoc(1,fltarr(512)) 


openr,2,’[walden]slit.dat’ 

sl=assoc(2,fltarr(11)) 

slit=sl(4) pec@invists wl to 2972 vice 2972.875 
hold=fltarr(512) 


hold=indgen(512)*0.0 
mold(375)=1.0 


@at(0)=convol(hold,slit)/3.133 ;spectra=l/wl bin so that scale factor is 
;absolute intensity 


close,2 
meLose,1 


return 
end 
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s;Walden,1991 
pro Olt eZ4 vealc 


openu, 1, GlleZdi02 aa 
dat=assoc(l fltareiol. 


openr,2, |[walden|sitatadat] 
sl=assoc(2,fltarr(1l1)) 
slit=sl(1)} 
hold=Eltarr(512) 


hold=indgen(512)*0.0 
hold(214)=1.0 


dat(0)=Cenvoliiholarc ieee 
close,2 
close; i 


return 
end 


;program calculates synthetic spectrum 
;for O+ atomic emission @ 2470A 


s;this makes wl centered @ 2968.5A 


;spectra is 1/wl bin so scale faceeer 
fis absolute intensity 


oe 


*Walden,1991 
;This program creates a syntheic model for the O++ atomic emission @ 2853A. 
meomwortt 2853calc 


Seenuyt, OLTI 2853.dat’ 
Mateassoc(l1,fltarr(512)) 


openr,2,’[walden]slit.dat’ 
Seo soc(Z,fltarr(1l)) 
slit=sl(1) 


mold=fltarr(512) 


hold=indgen(512)*0.0 
mold(337)=1.0 


Spike=convol(hold,slit) 

Bat(0)=spike/total(spike)/3.133 ;dat(0) is regular convolved spike 

spike(338:510)=spike(339:511) ;dat(1) is spike with adjustment for apparent 
~eosmepped pont in data: 

dat(1)=(spike/total(spike))/3.133 


gelose,2 
close,l 


return 
end 


oS 


eWalden,1g991 
;This program calculates a synthetic spectrum for the NII doublet 
;with the ratio based on the findings of Bucsela and Sharp (1989) 


pro NII_ 


close,2 
close,l 
return 
end 


doubcalc 
openu,1,’NII_ 2143doub.dat’ 
nii=zassoc(1,fltarr(512) ) 


openr,2,’[waldenJslit.dat’ 
Sl=assoc(2,fltarr(11)) 


hold=fltarr(512) 
hold tos) =-s6 72139.68A line 
intl=convel(hold,sl(1))) jw) 1s8in-cenbemror sea 


hold2=fltarr(512) 
hold2(110)=1.0 221.4:3 7 S5aiae 
intZ=convol(hold2,sl1(3));wl is 1/4 to the right of Center (cream ee 


int=intl+int2 

mii(l)=(int/total(imt) ) secs 
;total of spectral area is 1/wl bin, so 
;scale factor will be the absolute intensity 
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sWalden,1991 

;This program creates a synthetic spectrum of atomic emissions observed in 
*the March 1990 MUSTANG data. 

pro otherspec 


openr,1,’[waldenjJslit.dat’ 
slit=assoc(1,fltarr(11)) 


meonu,2,’otherspec.dat’ 
Spect=assoc(2,fltarr(512) ) 


spec=fltarr(512) 

numlines=20 snumber of contributing atomics 
;wavelengths determined from Bureau of Standards 
;listing of atomic emissions and NASA listing of atomic 
-emissions 





; Wavelength Species Pixel 
; mSs46 .0 NIII eS 
; ye5 2.6 OI! P23 
; ZAi0. 6 Cit 125 
; 2, 3 Sa ara) ARII BZ 
; 2384.6 OIV 187 
; 2445.6 Orn 207 
; 2490.3 NII 2 all 
; 2 UCI NII Ji Jad} 
; Beo>.8 O11 277 
; i098 NII IS, AL 
; 2741.3 OO Sid) dk 
; Zoe. 3 NII 347 
3 m7 3.6,29/74.6, NII multiplet 375-377 
f Pew o.7, 2915.9, 

; menor 9 ,2977.3, 

: moog. 4,2978.6 


peer ot 25,1/72,187,207,221,227,277,291,301,347,$ 
eee, 5/760,376,376,376,377,377) 


Bac=(10,2,1, 
3, 


Deano a oO, 4, 4,2,3,5,4,4 
mno=(3,4, ’ 


H! 
Sewage, 1 ,1,1,1,3,2,4,1,3 


~ os 


2 

ror i=0,numlines-1 do begin 
blank=fltarr(512) 
Peeank(pix(i))=.1*int(i) 
impulse=convol(blank,slit(sino(i))) 

| spec=spect+impulse 

2ndfor 

| 

spect(0)=indgen(512)*3.133+1798 

spect(1)=(spec/total (spec) )/3.133 "Southat scale vtactor 1s absolute int. 

slose,2 

mose,l 


-eturn 
na 


5) 


APPENDIX C 


Appendix C is the programs used to determine the standard 
deviations of the means(sdoms) for the data. The sdoms are 
used to determine the weighting factors for the fitting 
routine by the program titled "weights2," which is also 


included. 
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swalden,1991 
oro m90sdev,alt,indexl,index2 


sThis procedure calculates the std.dev and std dev of the mean 

»for the averaged spectra. 

memeieencOnsist Of the requested altitude (1-25) and its corresponding 
sbeginning and ending index numbers. 

sOutputs are the SDOM (unshifted and shifted) and are written to 

;um90 sdom.dat and m90_ sdom.dat, respectively 


1avg=0 

st=// 
sigma=dblarr(512) 
sigmam=dblarr(512) 


gpenr,1,‘dub0:(claytonJmustang raw.dat’ 


Mat=assoc(l,intarr(512)) ;Opens the raw data file 
ypenr,2,’[mustang]um90 avg.dat’ ;opens the averaged spectrum file 
avg=assoc(2,fltarr(512) ) ; (unshifted version) 
mpenu,3,’{mustangJum90 sdm.dat’ ;opens the file to write SDOM into 


Sdom=assoc(3,fltarr(512)) 


1=index2-indexl 
For j=0,n do begin 
print, indexli+j 
Pant,’ DO you want to include this record’ 
meant,’ Type <CR> for yes, N for no.’ 
read,st 
if st eq ‘’n’ then goto,labell 
Sigma=sigma + (dat(indexl+j)-avg(alt))~2 
navg=navg+l 
labelli:dummy=0 






sndfor 

‘ree=navg-1l 
slgma=Ssqrt(sigma/free) 
slgmam=sigma/sqrt(navg) 
‘dom(alt)=sigmam 
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;Walden,1991 
pro weights2 


;This program checks for values where the detector saturated, and gives 
;weights based on the deviations of the mean being equal to the data at 

;those points. It also gives a small weight factor to those elements 
;corresponding to nitrogen second positive wavelengths (sets sdom equal to 
;largest data value in N2 2nd pos bands). Since sdom was given value of data at 
;dropped points, it gives small weights to dropped values (with the exception 
;0f 2972A) 


openr,l,’[mustang]m90 intsdm.dat’ 
dom=assoc(1,fltarr(512) ) 


openr,2,’([mustang]m90 int.dat’ 
dat=assoewe,: | tcammgat )) 


openu,3,’'m90 domwt.dat’ 
domwt= Aeeoets, f£ltarrC> ieee) 


openr,4,’{mack]mack cal.dat’ 
ss=assoc(4,fltarr(512) ) 
sens=ss(22) 


openr,5,’[mustang]m90 avg.dat’ 
dd=assoc(5,fltarr(512)) 


big=1800 ;approx N2 2nd pos largest datapoint value 
holddom=fltarr(si72) 

holddomwt=fltarr(512) 

holddat=fltarr(512) 

sat=900.00/sens 


for i=1,25 do begin 
holddom=dom(1i) 
holddat=dat(i) 
for j=0,511 do begin 
;if detector was saturated, make sdom equal to data value 
if (holddat(j) gt sat(j)) then holddom(j)=holddat(j) 
;if wavelength is in N2 2nd pos, make sdom equal to big wavelength 


if (j eq 277) or (j eq 278) or (j eq 279) then holddom(j)=big 72669 
if (j eq 318) or (3 eq 319) or (j) eg 320) then Nolcdon =529 ;27 96 
if (j eq 322) or (j eq 323) or () eg 324)) then holddom(y)—erg Zo. 
if (j eq 362) or (j eq 363) or (j eq 364) then holddom(j)=big 72934 
L£ (j eq 368) or (j eq 369) or (j eq 370) then holddom(j)=big 92936 
L£ (j eq 374) or (j eq 375) or (j eq 376) then holddom(j)=big AS) a) S 
L£ (j eq 409) or (j eq 410) or (j eg 411) then holddom(j)=big ; 306s 
L.£ (3 eq 417) or (3 eq 418) or (j eq 419) then holddom(j)=big sisi) 7) 
-£ (j eq 425) or (j eq 426) or (j eq 427) then holddom(j)=big »3132 
.£ (j eq 433) or (j eq 434) or {3 ©€q 435) then holddon( a eri, ree ales) 7: 
-£ (j eq 461) or (j eq 462) or (Jj eq 463) then holddom(j)=big 73244 
-£ (j eq 470) or (j eq 471) or (j eq 472) then holddom(j)=big 73204 
-£ (3 eq 480) or (j eq 481) or (j eq 482) then holddom(j)=big 7 330 
£ (j eq 490) or (j eq 491) or (j eq 492) then holddom(j)=big poe) /| 
£ (j eq 501) or (3 eq 502) of (3 Cg 503) then holddonya)—bagq P3308 
F (j gt 493) then holddom(j)=big ;last spike 


now give weights their values 
if holddom(j) eq 0. then holddomwt(j)=1. 
98 


1f holddom(j) ne 0. then holddomwt(j)=1./holddom(j)*2 
endfor 
sCorrect for where low value of data point at dropped word gives false 
;high weight to point by giving dropped points 1/10 the maximum weight 
;0f any other point 
meexcept for 2972A -pixels 371 to 376 and 2850A -333 to 336). 
dropped=where(dd(i) eq -1.) 
hold=holddomwt 
hold(dropped)=0.0 
holddomwt(dropped)=.1*max(hold) 
: holddomwt (371:376)=max(hold) 
: holddomwt (333:336)=max(hold) 
;assign weight to file m90 domwt.dat 
domwt(i)=holddomwt 
endfor 
close,l 
close,2 
close,3 
close,4 
close,5 


return 
end 


| 
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s;Walden,1991 

;This program shifts the MUSTANG data array according to dropped pixewcueae 
;input in the procedure shiftarray. 

pro shiftdata 


openr,1l,’{mustangJum90 avg.dat’ 


mm=assoc(1,fltarr(sl2)) 


openu, 2,” (mustang |mgUvavd aa. 
dd=assoc(Z2,fltarrt(o.2)) 


Opent, 3, [mustang | unvemedmedar” 
UU=aSsSoOc(S, fltarr( olZ)) 


openu,4,’[{mustang]m90 sdm.dat’ 
SS=assoc(4, fitaprr( 512) > 


wl=indgen(512)*3.133+1798 


dd(0)=wl 
ss(0Q)=wl 


for i=1,25 do begin 
shiftarray,mm(i),holdl 
dd(i)=holdl 
Shittartay subi) Nelda 
ss(i)=hold2 

endfor 


close,4 
close,3 
close,2 
close,l 


return 
2nd 
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;Walden,1991 

@uimemprogram shifts an array of data according to the "drop" array of 
;pixels. 

Meo shifttarray,inarr,outarr 


Setter 46, 430,414,398, 382, 366,350, 350,334,318,302,286,270,254,239,223,223,$ 
eee eel ores g, 145,127,111,95,95,80,80,64,48,32,16] 


ndrops=n_ elements(drop) 
outarr=inarr 
length=n elements(inarr) 


for i=0,ndrops-1 do begin 

; This section is to put a -1 in dropped values and shift as required. 
outarr(drop(i)+l)=outarr(drop(i):length-2) 
outarr(drop(i))=-1.0 

endfor 


return 
end 
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;Walden,1991 

;This program takes the MUSTANG data from the averaged, shifted records 
;and puts in the appropriate average, and adjusts for instrument response. 
pro average 


openr,1,’[mustang]m90 avg.dat’ 
mm=assoc(l,£&ltarr(sizy 


openr,2,’[mustang]m90 sdm.dat’ 


ss=assoc(2,ftltarr(s!zZ)) 


openu,3,’[mustang]m90 int.dat’ 
dat=eassoc(3,fltarr( 5122) ) 


openu,4,‘[mustang]m90 antsdmedat’ 
dev=assoc(4,fltarr(512) ) 


openr,5,’[mack]mack cal.dat’ 
sen=assoc(5,fltarr(512) ) 
sens=sen(22) 

boolean=fltarr(512) 

dat(0)=mm(0) 

dev(0)=ss(0) 


for i=1,Z5edeebeqin 
;get array of where dropped words are 
points=mm(i) 
index=where(points eq -1.0) 
;correct for sensitivity 
points=points/sens 
;assign average for dropped data values 
points(index)=(points(index+1)+points(index-1))/2. 
;correct averages for double drops 
boolean(index)=1 
double=where(boolean*shift(boolean,-1) ) 
points(double)=(2*points(double-1)+points(double+2))/3. 
points(double+1)=(points(double-1)+2*points(double+2))/3. 
;assign value of data to dropped word sdom’s 
errors=ss(i) 
errors(index)=points (index) 
zero first six elements 
points(0:5)=0.0 
errors(0:5)=0.0 
assign data to mgUmint aa 
dat(i)=points 
correct for sensitivity and assign sdom to m90 intsdm.dat 
dev(i)=errors/sens % 
ndfor 


:lose,5 
‘lose,4 
lose, 3 
OSC AZ 
lose,l 


eturn 
nd 
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piataen, 1991 
pro units 


Bpenr,1,’[mustang]m90 avg.dat’ ;reads out avg spectra, shifted 
a=assoc(1,fltarr(512)) 


openr,4,’{mustang]m90 sdm.dat’ ;reads out SDOMs, shifted 
c=assoc(4,fltarr(512)) 

Oopenr,2,’{mack]mack cal.dat’ fool obatlOumeGrcOmhect UnleS omiuntensity 
Beassoc(2Z,ftlitarr(512)) 

sens=b(22) 

Soenw,3,’(mustang]m90 int.dat’,512*4 ;file to accept corrected values 


mit=assoc(3,fltarr(512)) 


openw,5,’[{mustang]m90 intsdm.dat’ ,512*4 


Bmeaom—assoc(5,fltarr(512)) 


BOr i=1,25 do begin 
miei) = a(i)/sens 
intdom(i)=c(i)/sens 
endfor 


close,l 
close,2 
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;Walden,1991 
:This program assigns zeros to the first ¥sax values "of an array oreoe 
PprO ZeCrosix, Ename 


openu,l,fname 
a=assoc(1,fitarr(512) ) 


temp=flitarrgotz) 


for i=1,25 do begin 


temp=a(i) 
temp(0:5)=0.0 
a(i)=temp 

endfor 

close,l 

return 

end 
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APPENDIX D 
Appendix D shows the programs used to convert the scale 
factors obtained from the fitting routine to absolute 
intensities. The equations were derived from Barth (1965). 
The equation used for converting the scale factor to 


intensity for the VK(0,5) band is: 
eee ALGO a homie CG), 


where Ip 5 is the band intensity in Rayleighs, VKCALC(0,5) is 
the total intensity of the synthetic VK(0,5) band in 
Rayleighs, and SF(alt) is the VK scale factor for the altitude 
desired. The factor VKCALC is just the area under the 
synthetic spectrum curve for the transition and is obtained 
using the IDL command TOTAL() and multiplying by wavelength 
bin size, as seen in the program VK_ABSINT. 

The conversion for the LBH scale factor included all 
bands, and is performed in the program LBH ABSINT and the 
subroutine INTENCALC. The basic equation for this progran, 
derived from equations (4) & (8) of Barth (1965) is given by: 


LBHCALC(J,k) Yo Qj yV5, 


aliv’ 
TT Si 
V5,K FWj,0Fj,k 


where I is the total intensity of LBH emissions in 


LBH 
Rayleighs, LBHCALC(j,k) is the area under the synthetic 


spectrum curve for a transition band (j,k), crn is the Franck- 
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Condon factor for the (j,k) band, and vy is the corresponding 


frequency for each transition. The summation over v" is to 
include all possible upper states. The result of this 
equation was essentially the same for any (j,k). This thesis 


utilized the (4,12) transition. 


O'S 


*Walden,1991 
pmeeeDroOgram converts the scale factor from the nimpfit all routine into 
san absolute intensity for a particular transition band of Vegard-Kaplan. 
*Stored in vk absint.dat, the values are for the (0,5) transition in record l, 
mememeo,6) in record 3, and the (1,8) in record 5. 
pro vk absint 
~ openr,10,’finalpara.dat’ 

p=assoc(10,fltarr(25)) 

temp=p(0) 

alt=p(21) 

vksf=p(13) 

delvksf=p(14) 


epenu,ll,’vk absint.dat’ 
g@re=ascsoc(l1,fltarr(25)) 
dat(0)=alt 


wl=indgen(512)*3.133+1798 
meld=fltarr(25) 
mernold=fltarr(25) 


for run=2,6,+2 do begin 
1f (run eq 2) then begin ;F-C factors are not currently used 
Le ;in the determination of intensity, 
a — 5) ‘Due are included for.aruture use 
: gl=5.9e-4 DG ieveieonen mene =i), (ly, 
; q2=1.69e-1 F-Gatactor Lon 0, 5) 
P26 
endif 
if (run eg 4) then begin 
vVi=( 
= 6 
; gl=5.9e-4 Mwo-Catactor tor (070) 
: q2=1.89e-1 "PC Lactome Loum, 6) 
tewik—2 0 
endif 
1f (run eg 6) then begin 
el 
vo 
ql=3.32e-3 eFC. factorezonu( 1,0 
: q2=7.87e-2 “H=G factor Lor (Alo ) 
twk=1.5 
endif 
for a=0,24 do begin 
int=indgen(512)*0.0 
vksyn,temp(a),vl,v2,wl,int 
intz=int*twk 
held tay=—toraliimt )~3.155*vkst(a) 
errhold(a)=total(int)*3.133*delvksf(a) 





endfor 

dat(run-1)=hold 
| dat(run)=errhold 
: endfor 
Bf lose,10 

lose,1l 

earn 

nd 


[synspec]vksyn 
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>Walden,1991 
pro Tbhmeboiae 
; This program determines the absolute intensity of LBH emissions based 
son the fits of intensity to three different transition bands (the results were 
;basically equal in Walden’s analysis). Stored in record 1 of lbh absint.dat 
sis based on th (2,9), record 3 1¥S8" the (4,12)— recorsa.5 1s 5) 0s 
openr,10,’finalpara.dat’ 
peassoce(l0,fltaruw~zs)) 
temp=p(0) 
lbhsf=p(9) 
dellbhsf=p(10) 
alt=p(21) 


opent 1 Py bnetous aac. 
t=assoc( li, fltarrece. 
WMovsneenc sl, gua ((0)) 


openu,12,’'1lbh absint.dat’ 
i=assoc(12,fltarr(25)) 
i(0)=alt 


holid=fitarr (255 
@rrhnold=f£ltarr is) 


£Or runs io sedempeain 

if (run eq 1)8ehen begin 

lee 

v2=9 

endif 

if (run eg 2) then begin 

ee al 

a — ee 

endif 

if (run eq 3) then begin 

vis=3 

Vien 

endif 

for a=0,24 do begin 
intencalc,temp(a),vl,v2,bandint 
intlbh=bandint*3.133*lbhsf(a)/lbhtotint((temp(a)/100)-1) 
errintlbh=bandint*3.133*dellbhsf(a)/lbhtotint((temp(a)/100)-1) 
hold(a)=intlbh 
errheld(a)=errint len 

endfor 


i1(2*run-1)=hold ;put absolute intensities in appropriate record 
1(2*run)=errhold ;PUL ELF EOKS In coLnecEsbecord 
evoleraie 


-lose,10 
-lose,ll 
lose,12 


-eturn 


nd 
[synspec]lbhsyn 
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»>Walden,1991 

mummomerocedure is called by the lbh absint program in order to determine 
mememmmtensity under a particular "hump" (a transition band) in the LBH curve. 
-The inputs vl and v2 are the upper and lower states’ vibrational levels 
;for the transition to be determined. Intlbh returns the intensity value. 
meemrmeencalc,T,vi1,v2,intlbh 


openr,l1,’[synspec]frnkcondn.dat’ 
fcxassoc(1,fltarr(15,30)) 
Pohic=fc(0) 


Seenr,2, (|walden]lbhtrans_ wl.dat’ 
Peea=assoc(2,fltarr(15,30)) 
lbhwl=wl(0) 


summa=0.0 
wl=indgen(512)*3.133+1798 
mye=tcitarr(512) 

Mont req=titarr(15, 30) 


for j=0,14 do begin 
for k=0,29 do begin 
mao lOnwl(i,k) ne 0) then lbhfreq(j,k)=2.997e8/lbhwl(j,k) 
mem tohwl(j,k) eq 0) then lbhfreq(j,k)=0.0 
endfor 
endfor 


mor 1-0,29 do begin 
summa=summa+(lbhfreq(vl,i)°3)*lbhfc(vl,i) 
endfor 


measyn,l,Vi,v2,wl,int 

mntviv2=total(int) 

eae lon=intviv2*summa/(lbhfreq(vl,v2)°3*lbhfc(v1,0)*lbhfc(vl,v2)) 
close,l 

close,2 

return 


end 
@{synspec)lbhsyn 
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APPENDIX E 


Appendix E contains plots of the final fits obtained, 
along with the program used to print them. The plots are 
divided into 400A increments, from 1800A to 3400A, and are in 
order of increasing altitude. Each plot shows the fit versus 
data in the upper view, and a breakdown of the fit by major 


contributors in the second view. 


Jia, 


»>Walden,1991 
pro plotapp,alt,xmin,xmax 


Seent,l,’ tinalpara.dat’ 
para=assoc(1,fltarr(25)) 


height=para(21) 


openr,2,’{mustang]m90 intnh.dat’ 


meeeecoc(2,fitarr(S512)) 


Siemonow Cinaltwk.dat’ 
Bereeagews,;tltarr(51i2)) 


temp=para(0) 
it=(temp(alt-1))/100 


Seene,4,’allvk.dat’ 
me—asesoc(4,fltarr(512)) 
pvk=para(13) 


Seemrno, all iIbh.dat’ 
meneeassoc(5,fltarr(512)) 
plbh=para(9) 


ee@emr,o, nodlt.dat’ 
See aesoc(6,fltarr(512)) 


Seentyy, noeps.dat’ 
Seem es5oc(/,fltarr(512)) 


openr,8,’nogammack.dat’ 
gam=assoc(8,fltarr(512)) 


opno=para(1ll) 


Meemt,7,’OIIL 2853.dat’ 
mmmv—assoc(9,fltarr(512)) 
poiii=para(7) 


Meont,19,’OII 2470.dat’ 
Smp—assoc(!10,fltarr(5l2)) 
poii=para(5) 


Meenmr,li,’OL 2972.dat’ 
Mr=acsoc(li,fltarr(512)) 
poi=para(3) 


meeticr,i2,’NII 2143doub.dat’ 
Me=assoc(12,fltarr(512)) 
pnZ=para(is) 


wl=indgen(512)*3.133+1798.000 


WeeodEanveplOes Gata VS. fLiltvand fit by 

;parts in format used by Walden thesis 

TEGimap pendix alt Canebes!l to 25, xmin 
;and xmax are presumed to be 1800-2200, 
pe2UWeZOCO WZ60G—=30007 or 3000-3400, so 

;Other combinations would need adjust- 

;ments to the plot key 


data 1S ln recoras G=lbato, 25 


JLieSeace IM Eecords f=#l to 25 


reel 


Set plotpanp: 

device, /portrait,/plotter on _off 
device, xof&tsct—-lac7 

device, yoffset=.635 
device,xsize=17.78 
device,ysize=24.13 


;TOP OF PAGE--DATA AND FIT 
set xy,xmin,xmax 
S@C VieGWDO tye, boy cc eee 


if (alt ge 18) then $ 

Imtitle=’Data & Fit at ’+strtrim((fix(height(alt-1))),2)+’km (downleg)’ 
if (alt. Ft _ 18) then. > 

tmtitle=’Data & Fit at ‘’+strtrim((fix(height(alt-1)))>2)2’km (upleque 


txtitle=’Wavelength(A) ’ 
'ytitle=’Intensity(R/A) ’ 


dat=d(alt) 
Elt=—f Cae 


!linetype=0 
plot, wiljaar 


tlinetype=2 
Op UGiey wal yemaste 


;KEY FOR TOP PLOT--xrange must be determined in terms of pixels so that 
;yrange can be known for placing the key in the correct place on the graph 
;case section does this for the ranges designated 


case l of 

(xmin eq 1800) and (xmax eq 2200): begin 
ldat=dat(0:129) 
PEL t=fre 0-17) 
end 

(xmin eq 2200) and (xmax eq 2600): begin 
Idat#dat( 1292257) 
Fiirt=f1t4 23 2 
end 

(xmin eq 2600) and (xmax eq 3000): begin 
ldat=dat( 257-307) 
Lilitefret 257-364) 
end 

(xmin eq 3000) and (xmax eq 3400): begin 
ldat=dat(384:511) 
Dire=f£i tie. > i) 


end 

else: begin ;key will probably not be correct 
ldat=dat ;this section must be manually altered 
lfit=fit ;to concur with the wl range 
end 


endcase 


;this section designates x and y coordinates for the sample lines in the key 
xcoordl={ xmin+20 .xmin+ 0) 

xcoord2=(xmin+20,xmin+100] 

xcoord3=(xmin+120,xmin+200] 

ycoord1l=( .96*max(ldat),.96*max(ldat) } 
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A aaa es ae 90*max(ldat) 
megeras—(|.84*max(ldat),.84*max(ldat) 
meoord4=(./7/8*max(ldat),./78*max(ldat) 
megords=(,./2*max(ldat),.72*max(ldat) 
memermdo=(|.66*max(ldat),.66*max(ldat) 
meoord7—(|.60*max(ldat),.60*max(ldat) 


!linetype=0 

Sect, scoordl,ycoordl 

tlinetype=2 

moor, xcoord!,ycoord2 

t!linetype=0 
mpoues,xmin+60,.995*ycoordl(0),’Data’ 
meomiecwxmint60, .995*ycoord2(0),’Fit’ 


SeeeeormeOr PAGE--PLOT OF PARTS AND FIT 


set xy,xmin,xmax 
Meee wpOrt,.18,.82,.09,.41 


mimctitle=’Data and Fit by Parts’ 
Ixtitle=’Wavelength(A)’ 
fytitle=’Intensity(R/A)’ 





f!linetype=0 
mbou,wl,dat 
moplot,wl,dat pmakes Gata a bold liane 
t!linetype=l1 
moplot,wl,vk(it)*pvk(alt-1l) ;dotted 
tlinetype=2 
Bevoe,wl,lbh(it)*plbh(alt-1l) ;dashed 
f!linetype=3 
memerbe=-pno(alt—-1)*(gam(it)+eps(it)+.25*dlt(it))/3.133e6 
mploc,wl,nitric >dash-dot 
tlinetype=4 
Bebmot,wl,0111(0)*poiii(alt-—l) ;solid 
tlinetype=5 
oplot,wl,n2(1)*pn2(alt-1) ;solid 
tlinetype=5 
oplot,wl,o011(0)*poii(alt-l) Acie Jeatie| 
tlinetype=5 





oplot,wl,oi(0)*poi(alt-1) SOG 
: 


key for bottom plot--x and y ranges from top plot used for locating key 


case 1 of 
mexemin eq 1800) and (xmax eq 2200): begin 
flinetype=0 
moon, xcOord2,ycoordl 
Selon, xcoord2,ycoordl 
!linetype=3 
| @pet, xCOOLd2,ycoord2 
flinetype=2 
| DmLow, XCOOLGL,VCoord 3 
| f!linetype=5 
Roloc  KeOOrd2, ycoord4 
!linetype=0 
mZoOMeo cnn lO IF o*yocoord!(0),’Data’ 
a Oueuc, mint nO 7. 99 S*veoqord2(0),‘NO’ 


————————w 


xyouts,Xmin+110, .995*veoerds( 0) eee 
xyouts, xmin+l10,.995*ycoord4(0),’Nit deus tome ee 
end 
(xmin eq 2200) and (xmax eq 2600): begin 
!linetype=0 
ODO, xCoenrd., Vooora: 
Oplot, xcoonadZ ; Veco ua 
flinetype=3 
oplot, xcoeerd2Z ,yveeoudZ 
f!linetype=2 
oplot, xcoord2, Veoouas 
!linetype=1 
oplot,xcoord2, ycoouds 
f!linetype=5 
oplot, xcooLrdZ, ycoomea. 
!linetype=0 
xVouts,xman+l 107. 99S*Vcoordi Ut Pac ae 
XVOuts ;-xmin+l10, .995*ycoordz(0) 7 23 
xyouts, xmin+l10,.995*Vcecrds (0) rare 
xyouts, xmin+ll0;7 .995*Vcoordst. Vi 
xVOuts,kxmintll0,.995*Vycocordsa( 0) Otte ee 
end 
(xmin eq 2600) and (xmax eq 3000): begin 
!linetype=0 
Oplol, xcoords,yeeocd: 
OplOot,xcoords, yeooral 
f!linetype=3 
oplot,xeoords ,vecord. 
flinetype=l 
ODLOL, MeOOLds, Woe mar 
f!linetype=5 
ODLOt, xcOOrds, yecooncac 
!linetype=4 
oplot,xcoord3s, yecoords 
flinetype=0 


xyouts,xmin+210, .995*ycoordl ‘Data’ 


(0), 

XVOUtTS , xmin+2Z10, 5995" 7ecoerd0 jae: 
xvyouts,xmin+210 7. 995*veoerds aie VK" 
XVOUtS, xMin+2109e 9IS*Vcoord4 (0 7 Oley 2A 
XVOUtCS,xmin+Z10, ,995*ycoords (0) Orla cose. 
end 

(xmin eq 3000) and (xmax eq 3400): begin 
!linetype=0 


Oplet, “cOOlaZ, ,eoord! 
Oplot, xcOOrd2, yeoonral 
!linetype=3 
Oploty, *GeordzZ yecordc 
flinetype=l 
Oploe, xcCOgnrd. 7 ceoras 
f!linetype=0 


XVOuts, xmint+110, .995*vyecoordi( 0] vaca. 
xyouts,;xmin+l10, -995*vcoerdZ( 0) Te 
xyouts,xmin+110). 99oeycoords Gin. 
end 

2lse: dummy=0 


=ndcase 


-close, 13 
close,12 
close,ll 
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elose,10 


 close,9 


close,8 


 close,?7 


close,6 
close,5 
close,4 
close,3 
close,2 
elose,l 


Return 
end 
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